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SUMMARY 


This work contributes to the study by the National Aeronautics and Space 
Administration (NASA) of the possibility that precipitation scavenging of HCi 
from large solid rocket (SRM) exhaust clouds may lead to unacceptably acidic 
rain in the Cape Canaveral, Florida, area before atmospheric dispersion has 
reduced the potential hazard to safe limits. Washout coefficients are derived, 
a previously developed HCil(g) washout model for rain scavenging of SRM clouds 
is refined and generalized, and the model is applied to nine independently 
determined SRM cloud dispersion cases. 

Derived washout coefficients A = AH*^, given in sec"^, are summarized in 
terms of rainfall rate H, given in mm rain/hr, and empirical constants (A,b) 
as follows: Washout of HCil(g) was based on the modified Frossling equation for 
convective diffusion of HCil(g) to a falling droplet; (A = 1.80 x 10“^, 
b = 0.565) applied for the Marshall-Palmer (M-P) raindrop size distribution, 
and (1.08 x 10~^, 0.625) applied for Kelkar's extensive raindrop size- intensity 
data for average rains. Washout of HCA (g + aq) was also derived, based on 
recent laboratory data on rain scavenging of diluted and humidified SRM exhaust 
aerosol; (1.52 x 10'^, 0.658) applied for the M-P distribution. Since this 
result stemmed from limited data and was statistically equivalent to A for 
HCil(g) washout, it is provisionally recommended that the geometric mean A for 
HCjl(g), i. e. (1.39 x 10“^, 0.595), be used to characterize HCA washout from SRM 
clouds. 

The washout model treats the idealized case of an independently generated 
vertical rainfall that overrides and scavenges an independently advecting and 
dispersing SRM cloud under stable stratification conditions in the lower tropo- 
sphere. The resultant pH of rain and HCil deposition rate are characterized in 
terms of (a) initial source strength of HCil; (b) dispersive decay of vertical 
HCil(g) column density a expressed in terms of case-specific empirical con- 
stants a and B and downwind distance (time), vis, a = aX”^; and (c) mean 
wind speed, downwind location of rainfall onset. A, and H. Cumulative areal 
deposition of HCA is characterized in terms of ground coordinates when 

(a) to (c) are specified and the SRM cloud geometry is assumed to be an expand- 
ing right circular cylinder. 

The washout model was applied to a refined "spring fair weather" (SFW) 
Space Shuttle case and eight Titan III (60 percent less exhaust) dispersion 
cases. These O decays were previously deduced by application of a multilayer 
Gaussian diffusion model to seven standard meteorological regimes for overland 
advection at Cape Canaveral. The Titan III dispersive decays of a and hence 
potential fffl, defined as volumetric average pi at first onset of rainfall, 
differed greatly among the seven regimes. A range of more than 2 pH units was 
spanned at X ^ 100 km downwind and t S 2 hr. At shorter distances the total 
span was less, but still exceeded 1 pH unit for X > 1 0 km and t > 0. 2 hr. 
Environmentally significant ffl's (^1.5) for infrequent exposures were shown 
possible at X S 50 km and t ^ 5 hr for the two most severe, least disper- 



sive Titan III cases. This result contrasts with a measured volumetric average 
pH of 4.61 ± 0.22 (monthly standard deviation) for rains at Cape Canaveral over 
the last 2 years. 

Detailed examples of downwind rainwater pH and for both potential and 

progressive washout, are shown for the least and most dispersive Titan III cases 
and the refined SFW Shuttle case. Both pH and increased with H, as a con- 

sequence of b < 1.0. High H occurring close to the launch site resulted in 
much greater G^ than lower H occurring farther out. Nearly ccmiplete HCA 
washout could occur within 30 km at high H. Although potential pB, applies 
strictly at rainfall onset, progressive washout at low H over large X 
resulted in slew pH change. Thus, damage to ground-receiver surfaces may tend 
to be greater (lower pH) and more extensive (longer footprint) at lower H. 

Several factors affect the validity of the predicted pB* s and HC?^ deposi- 
tions. First, the model application specifically excludes convective activity 
and, particularly, dynamics of rainout processes when HCA and alumina interact 
with natural clouds. This exclusion reduces the utility of the model, since 
convective activity occurs frequently at Cape Canaveral. Clearly a more real- 
istic atmospheric dynamics model is needed along with parameterizations of the 
essential cloud microphysics- scavenging processes. Second, uncertainties 
related to the dispersion calculations, and hence O decay inputs, affect the 
validity of predictions, especially at large X. Inclusion of (a) convective 
loss of HCil from the SRM cloud's upper boundary, (b) HCJl sorption at ground 
level, and (c) a more realistic treatment of horizontal wind shear effects 
would have increased potential rain pH. However, the values of o were based 
on variances derived from relatively small-scale turbulence measurements; vari- 
ances scaled more appropriately to cloud size would have resulted in smaller 
but more realistic values of 3 at large X and t, similar to those observed 
experimentally for recent Titan III clouds. This modification would decrease 
potential rain pH. Such competing effects require complex analyses, and their 
net effects are presently not predictable. Third, the model applies for ideal- 
ized HCA washout by rain of specified intensity and average polydispersity. 
Naturally occurring rains vary widely in both respects, and their temporal- 
spatial characteristics are only roughly predictable. Finally, the washout of 
HCil (g + aq) still remains an uncertainty. Although HC!^(g) will predominate over 
HCil (aq) on chlorided- alumina nuclei after significant cloud dilution at low- 
to-moderate humidities, the present values of A for HCil (g + aq) may be signif- 
icantly in error at very high humidities (approximately ^95 percent), which fre- 
quently exist at Cape Canaveral. 


INTRODUCTION 

NASA has been examining the possible environmental impacts of its Space 
Shuttle Program for about 8 years. Formal Environmental Impact Statements were 
published in July 1972 (ref. 1) and April 1978 (ref. 2). One of the cited 
potential problem areas, atmospheric pollution, stems mainly frexn use of a 
solid- propellant rocket motor (SRM) booster design. The tropospheric portion 
of this environmental problan centers on the possible effects of relatively 
large, localized, low-level releases of SRM exhaust products, which include 
more than 70 metric tons of hydrogen chloride (HCil) and about 110 metric tons 
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of aluminum oxide particles (0*01 to 20 ]jm) emitted per launch below 4-km alti- 
tude. More specifically, the possibility exists that precipitation scavenging 
of HCA may lead to localized deposition of unacceptably acidic (hydrochloric 
acid) rain on nearby land areas or protected waters before atmospheric disper- 
sion has effectively reduced the potential hazard to safe limits. 

One purpose of this paper is to refine and generalize a previously devel- 
oped model (ref. 3) which characterized the idealized washout of gaseous HCd 
and resultant deposition of acidic rain from an advecting and dispersing SRM 
cloud. A second purpose is to characterize the possible effects that various 
meteorological conditions at Cape Canaveral, Florida, may have on acidic rain 
deposition. 

Obviously the development of understanding and quantitative methodology 
for predicting HCil scavenging and for assessing possible environmental conse- 
quences of various acid rain deposition is an essential precursor to specifi- 
cation of safe, realistic launch constraints. The original Environmental 
Impact Statement (EIS) for the Space Shuttle Program (July 1972) states, 

”. . . for the (predicted) overland trajectories of the exhaust cloud, the 
possible harmful effects of rain containing HCA will be analyzed prior to each 
firing. If the calculations predict unfavorable conditions, the launch will 
be postponed." The recent revision of the EIS (Final, April 1978) states, 

". . . the results of this program ... to model the occurrence of acidic 
rain . . . and to predict the acidity of rainfall that might occur . . . will 
provide a model to define in advance the go-no-go (launch constraints) cri- 
teria to minimize unacceptable environmental effects from acidic rainfall." 

The absorption of HCil(g) by free-falling aqueous droplets was previously 
studied experimentally (ref. 3) . The data were shown to be well characterized 
the modified Prossling equation with an appropriate binary diffusion coef- 
ficient^ for HC?/ in air. Also, an idealized rain scavenging model was developed 
and used to characterize the washout of HCJl(g) frcxn SRM exhaust clouds (ref. 3). 
This model treated the simplified case of an independently generated vertical 
rainfall that overrides and scavenges by washout processes an independently 
dispersing SRM exhaust cloud. The model was developed to apply strictly at low- 
to-moderate SRM cloud relative humidities, where HCJl(g) tends to predominate 
over the aqueous acid aerosol ccmponent after a few minutes of cloud dilution, 
and under stable stratification conditions in the lower troposphere. 

The model in reference 3 was used to predict the resultant characteristics 
of acid rain deposition for an SFW Space Shuttle SRM exhaust cloud dispersion 


^Although the initially assumed binary diffusion coefficient of 

0.20 cm^/sec for HCi in N 2 (air) was shown to fit the four column-length — 
residence- time sets of droplet absorption data well and thus was used through- 
out the reference 3 analysis, the four data sets were best fit statistically 
^AB “ 0.175 cmVsec. It was subsequently found that use of a theoretical 
equation based on modern kinetic theory and the Lennard-Jones expression for 
intermolecular forces (ref. 4) led to a calculated value of 0.170 cmVsec at 
25^C and 1 atm. This latter value was adopted for the present paper and then 
corrected to 0.187 cm^/sec at 15^C and 0.85 atm. 
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case, derived independently from application of a Gaussian multilayer diffusion 
model, MDM-4(II), that was developed at the NASA Marshall Space Flight Center 
(MSFC) . The SFW dispersion case was based on one of seven standard meteorol- 
ogies documented for the Cape Canaveral area* These meteorologies were 
originally selected to represent a range of characteristic cases of overland 
advection with turbulent diffusion in the planetary boundary layer in order to 
establish an initial basis for assessing postlaunch HCil concentration histories 
at the Earth's surface (ref. 5). 

In the present paper, the idealized washout model is first refined and gen- 
eralized and then is applied to nine independently determined exhaust cloud dis- 
persion cases (ref. 6 and an unpublished paper by G. L. Pellett) which were 
derived from the set of seven standard meteorological conditions. The model 
improvements and additions consist of (a) including more representative values 
for air properties and HC^/ diffusion coefficient; (b) developing a more compre- 
hensive analytic approach for generalizing HCii^(g) absorption and washout and 
thus resultant rain pH and HCZ deposition characteristics; (c) deriving two new 
HCil(g) washout coefficient expressions based on the modified Frossling equation 
and integrated over both the well-known Marshall-Palmer (M-P) raindrop size dis- 
tribution and the previously used Kelkar raindrop size- intensity data; and 
(d) deriving a new HCi. (g + aq) washout coefficient expression for SRM exhaust 
aerosol based on a detailed review and analysis of published laboratory data 
on rain scavenging of solid rocket exhaust and integrated over the M-P distri- 
bution. While the previous washout model was applied only to the Shuttle SFW 
meteorological case, the refined washout model is applied to (a) an improved 
version of the Shuttle case, (b) all seven standard meteorological cases for 
the smaller but chemically similar Titan III SRM propelled vehicle, which 
exhausts «40 percent of that for the Shuttle in the altitude range 0 to 4 km, 
and (c) an abnormal (pad abort) Titan III case. The Titan III cases provide 
a basis for detailed comparisons with existing in-cloud data on Titan III 
launch effluents and also associated acidic rainfall in one case. Finally, 
the characteristics of acidic rainfall which might result frcan either vehicle 
under various meteorological conditions are clearly demonstrated; i.e., 
explicit analytic predictions of downwind rainwater pH and HCH ground deposi- 
tion are given for various meteorologies, assumed onsets of rainfall, and rain- 
fall intensities. 

Use of trade names or names of manufacturers in this paper does not consti- 
tute an official endorsement, expressed or implied, of such products or manu- 
facturers by NASA. 


SYMBOLS 

A empirical constant in power-law expression for washout coefficient 

(eq. (41)), sec"^ 

b empirical exponent in power-law expression for washout coefficient 

(eq. (41 ) ) , dimensionless 
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t 

C(HCil) concentration of HCil (g + aq) in diluted SRM exhaust, as used by Fenton 
and Purcell (ref. 25) and expressed in dimensionless volume/volume 
units, ppnv 

/s 

C(HCA) peak value of measured in-cloud HCiZ/ concentration, expressed in dimen- 
sionless volume/ volume units, pjrniv 

Cf molar air density, 3.60 x mol air/cm^ at 15^C and 0.85 atm 

^AB binary diffusion coefficient for HCA in N 2 or air; 0.187 cmVsec is 

used at 15^C and 0.85 atm (see footnote 1) 

Dx effective diameter of cylindrical (right circular) SRM exhaust cloud 

in horizontal plane, m 

d raindrop diameter, cm 

E(a,d) collection efficiency for falling drop of diameter d which collects 
aerosol particles and droplets of radius a; efficiency represents 
fractional collection of particles and/or droplets within geometric 
cross-sectional area TTd^, dimensionless 

Gx cumulative deposition of aqueous HCil (in rainwater) on ground, at 

location X along SRM cloud- centroid path, g 

Gx,y cumulative deposition of aqueous (in rainwater) on ground, at 

location (X,Y) along SRM cloud path, g 

Gx rate of deposition of aqueous (in rainwater) on ground, at loca- 

tion X along SRM cloud- centroid path, g HCA/hr-m^ 

H rainfall accumulation rate, mm rain/hr 

Hjjj upper bound of surface mixing layer, m 

kx gas-i^ase mass- transfer coefficient, mol HCil/cm^-sec 

M;^ mass concentration of HC?.(g + aq) in diluted SRM exhaust, as used by 

Fenton and Purcell (ref. 25) , g HCJl/m^ 

m(HC5/) molarity of aqueous HC5/ in droplets; mixing-cup average molarity for 
polydisperse rain, mol HCi/h 

niQ initial mass of HCil in SRM exhaust cloud, g HCS/ 

n number of raindrops per unit volume of air, an“*^ 

An^ incremental number of raindrops per unit volume for size class i, 

cm"^ 

An^ raindrop intensity for raindrops of size class i, rtT^-sec^^ 
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p (HCA) 


i 


average concentration of HCil(g) in parcel of diluted SRM exhaust, 
expressed in dimensionless volume/ volume units, ppmv 

Pq(HCJI) initial value of p(HCil) before HCA(g) washout begins at t = 0, 
ppmv 

pH measure of acidity due to H*^(aq); defined for an idealized solution 

of fully dissociated HCil{aq) by log*|o [l/m(HCA) ] 

Qj effective heat release for solid rocket propellant combustion, cal/g 

propellant (1 cal = 4.184 J) 

rate of HCA scavenging by individual droplets, g HCil/sec 

Sh Sherwood number based on droplet diameter (eq. (3)), dimensionless 

t elapsed time after onset of ground deposition due to UCH washout, sec 

tg exposure time for droplet absorption during fall through parcel of 

gas containing HCA, sec 

tQ elapsed time corresponding to Xq, hr 

t* transit time for droplet to fall to Earth's surface from altitude 

z, sec 

Uq mean unidirectional wind speed which describes horizontal motion of 

SRM exhaust cloud, m/sec 

Vqq droplet terminal velocity during free fall through parcel of air, 

cm/sec 

molar flow of HCil(g) to droplet surface during HCil(g) absorption, 
mol HCiz/sec 

X unidirectional downwind distance frcan launch site, Xq + X^g, km 

Xq virtual source distance upwind from SRM cloud stabilization point at 

^cs “ ^ ^ 

X^g unidirectional downwind distance from SRM cloud stabilization point, 

km 

Xj. unidirectional downwind distance from launch site where onset of uni- 

form rain occurs (before arrival of SRM cloud) , km 

Xao mole fraction of HCA(g) in equilibrium at gas-droplet interface 

XAco mole fraction of HCA(g) in bulk gas 
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Y perpendicular distance from cloud- centroid path, km 

Zjn SRM cloud-centroid height, m 

z vertical distance for droplet fall through SRM exhaust cloud, m 

^top height of SRM cloud top above ground, m 

a, 3 constants which define power-law decay of HCil(g) column density a 

in dispersing SRM exhaust cloud; a is equivalent to ppmv-m at 
X = 1 km, and 3 is dimensionless exponent of X in equation (48) 

A washout coefficient defined by equation (8) and generalized by equa- 

tion (41), sec”^ 

Vf kinematic viscosity of air, 0.172 cmVsec at 15^C and 0.85 atm 

Pair mass density of air, 1041 g air/m^ at 15^C and 0.85 atm 

a vertical HCil(g) column density, ppmv-m 

Subscript: 

pot potential 

Abbreviations ; 

CC cloud centroid 

CFP cold front passage 

F-P Fenton and Purcell 

FFW fall fair weather 

FW, Pre-CF fair weather pre-cold front 

LLSB low-level sea breeze 

MDM multilayer diffusion model 

M-P Marshall and Palmer 

Post-CF post- cold front 

SB sea breeze 

SFW spring fair weather 

SRM solid rocket motor 
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DEVELOPMENT OF WASHOUT MODEL 


The first two sections of this development follow closely that given in 
reference 3, but the coefficients have been revaluated to conform to the revised 
HCil(g) diffusion coefficient and the newly defined ambient conditions at 15^C 
and 0.85 atm, which more closely represent average conditions aloft at Cape 
Canaveral. Subsequent equations depend on these coefficients. The next section 
on washout coefficient for monodisperse rain serves a dual, although not immed- 
iately obvious, purpose. First, it defines the washout process in terms of 
droplet characteristics, and this leads to a set of analytic expressions that 
can be used independently to characterize washout by droplets of specified size. 
Second, it serves to define expressions that are needed later in the derivation 
of washout by polydisperse rain, when raindrop intensity data (e.g., Kelkar) 
are used to evaluate washout coefficients at various rainfall intensities. 


HQil Absorption and Resultant pH of Falling Droplets 

Hydrogen chloride is a highly soluble gas, and within an SRM cloud, HC?/(g) 
transfer to a falling droplet may be considered to be gas-phase diffusion 
limited. For low mass-transfer rates and where gas-phase mass-transfer resis- 
tance controls, droplet absorption can be described by (ref. 7) 


-Wa = 


Xao - 


1 - X 


Ao 


( 1 ) 


Now will always be much smaller than unity in an SRM cloud. Further- 

more, it will be shown later that, for all practical cases of HCil(g) washout 
(i.e., for ^0.5-mm rain/hr falling through a dispersing SRM cloud), X^^q will 
always be much smaller than This is due to the high solubility of HCA 

and to the expected range of HCil(g) column densities and concentrations. How- 
ever, note that, for the potentially reversible case of in-cloud scavenging, 

^Ao ^Aoo will tend to converge as HCl equilibrium is approached in cloud 

aerosol droplets, since they are much smaller and also have substantially longer 
residence times. Thus 


Wa = kx’f<3^XA^ 


( 2 ) 


The Frossling semiempirical mass-transfer correlation (ref. 8) , as modified 
by Ranz and Marshall (ref. 9) , has been successfully applied to a number of 
cases of mass transport to (or from) a sphere immersed in a moving fluid 
(ref. 7) . The dimensionless Sherwood number (diameter basis) is given in terms 
of the Reynolds and Schmidt numbers as 

M /‘3V \ 

Sh = = 2 + 0.60 (3) 

°f°AB \Vf y \Dab/ 
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By taking values for the properties of air at 15°C and 0.85 atm (see symbols) 
and adopting a value for the binary diffusion coefficient (D^ = 0.170 cmVsec 
at 25®C and 1 atm, corrected to = 0.187 cmVsec at 15°C and 0.85 atm; see 

footnote 1) , equation (3) becomes 


Sh = 2 + 1 .407(dV„) V2 


(4) 


Recognizing that p(HCil) = 10®X^^ and substituting for and in 

equation (2) , a pair of expressions^ is obtained for the molar absorption rate: 


Wa = 10-6cfD^TTd(Sh)p(HCil) (5a) 

Wa = 10~^l(4.23d + 2.98d3/2vJ/2^p(HCil) (5b) 


The resulting molarity of HCil in a falling droplet after an absorption exposure 
time t^ is 


m(HC«-) = 6 X 10"3cfD;^3d-2(sh)p(HCA)te (6a) 

m(HCil) = 10"8(8.08d"2 + 5.68d-3/2vJ/2)p(HC«,) t^ (6b) 


Assuming that the absorbed HCil is 100 percent dissociated to H'*', ignoring 
for a moment the initial H*** concentration in a rain droplet, and taking the 
exposure time for a droplet falling at terminal velocity through a cloud of 
dilute HCA(g), z meters thick, te = lOOz/V^ (no net vertical air motion in 
the SRM cloud) , the final pH of a fallen droplet obtained is 


pH = log 


pH = log 


1 0d2v„ 

CO 

6CfDABP(HCS,)z(Sh) 

10®d2v 


(8.08 + 5.68(dVJ V23 p(HCA)z 


(7a) 

(7b) 


Applicability of Irreversible Absorption Assumption 

The applicability of equations (2) to (7) to an SRM cloud is now checked 
for a typical case. An analysis of the SFW Shuttle exhaust cloud disperson 
case, described later, led to predicted values of p(HCA) and z at various 


^The second expression is included to indicate d and dependencies 

and also to facilitate user canputations. 
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distances from the launch site. Two sets of values within the main region of 
interest, 10 to 100 km, have been used in equation (7b) to calculate pH and 
molarity for rain droplets in the size range of interest, 0,3 to 0.03 cm. 
Results are shown in table I along with values corresponding to equilib- 

rium vapor pressures of the respective HC^t solutions at 15^C and 0.85 atm 
(calculated from ref. 10 vapor-pressure formulas). 


TABLE I.- SAMPLE GALOJLATIC^JS AND CHECK ON APPLICABILITY OF EQUATIC^I (2) 
FOR DROPLETS FALLING THROUGH AN SRM EXHAUST CLOUD* 


d, 

cm 

Vco' 

cm/sec 

pH 

(100% ionized 
to H+) 

m(HCH) = 10“PH, 
mol HCA/L 
in droplet 

XAco» 
HC)l in 
bulk gas 

Xao/’*' 
HCS, at 
gas- liquid 
interface 


PQ(HCi.)z = 2 X 10^ ppmv-m at X = 10. 

9 km downwind distance* 

0.3 

860 

1.59 

0.0257 

5 X 10-6 

0.0001 X 10-6 

.1 

430 

.676 

.21 1 

5 

.0057 

.05 

218 

.006 

.986 

5 

.140 

.03 

123 

-.536 

3.44 

5 

7.94 


Pq(HCA)z = 10^ ppnv 

-m at X = 67.5 

km downwind 

distance* 

0.3 

860 

2.89 

0.0013 

0.25 X 10-6 

A 

o 

. 

o 

o 

o 

o 

X 

o 

o^ 

.1 

430 

1 .98 

.0105 

.25 

.00002 

.05 

218 

1.31 

.0493 

.25 

.00035 

.03 

123 

.76 

.172 

.25 

.0038 


♦Based on fully modified SFW Shuttle exhaust cloud dispersion case 
developed later in this paper: Po(HCil)z = lO^X”^*^"^ and mQ = 61 x 10^. 

’^X^o is based on empirical fit (ref. 10) of equilibrium HCS/ vapor- 
pressure data for 15^ aqueous ECl and 0.85 atm. 


In general, the calculated pH values in table I indicate acidities less 
than 1 mol HCil/L. The obvious exception is for droplets smaller than 0.05 cm 
falling at a close distance (11 km) to the launch site. Although the basic 
assumption used in deriving equation (2), i.e., X^oo ^Ao' is valid for 
0.05-cm droplets falling at 1 1 km downwind, it fails for droplets of smaller 
size (e.g., 0.03 cm) that fall close to the launch site. However, at larger 
distances (or alternately, with initially lower HCA column densities and/or 
shortly after washout begins) , the irreversibility assumption is adequate for 
droplets as small as 0.03 cm; at a distance of 68 km, the assumption is highly 
valid for all droplet sizes. 

Thus, for this range of raindrop size, taking into consideration the rela- 
tive contributions of various droplet sizes to integrated HCjl(g) washout (dis- 
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cussed later ) , the use of equation (2) in the present model is considered ade- 
quate for rainfall rates >0.5-iran rain/hr at all relevant distances (S10 Ion) from 
the launch site or, alternately, whenever HC& column density is smaller than 
2 X 104 ppmv-m. 


Washout Coefficient for Monodisperse Rain 

In order to describe the effective transient HCfc(g) concentration in an 
SRM cloud t seconds after the SRM cloud encounters a steady overriding rain, 
a solution to the scavenging- rate equation for irreversible washout (refs. 11 
to 14) is required. Thus 


dp (HCil) 
dt 


-Ap(HCA) 


(8) 


The washout coefficient A can be seen to satisfy its operational definition 
as the ratio of scavenging rate to the concentration of the scavenged component. 
The solution of equation (8) which applies to a finite thickness z (m) of 
the SRM cloud is 


p(HCil)z = pQ(HCil)z exp (-At) 


(9) 


For idealized monodisperse rain, of droplet size d and number density n, it 
can be shown from (a) the operational definition of A and (b) rearrangement 
of equation (5a) that the washout coefficient is simply 


A = = 7TDaRd(Sh)n (10) 

p(HCil)Cf 


Since the monodisperse rainfall rate is defined by 


H = eOOOTId^Vji 


( 11 ) 


rearrangement of equation (11) and substitution for n in equation (10) yields 
a pair of alternative expressions for the washout coefficient for monodisperse 
rain: 


HDAB(Sh) 

A = (12a) 

6000 d 2 v„ 

oo 
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(12b) 


A * 3.12 X 10“5H^2d“2v^~l + 1 . 407d“3/2v^“l/2^ 

Noter from equation (12a), that A is essentially proportional to since 

Sh is only weakly dependent on Dp^ in equation (3) . 

Since equations (7) and (12) contain (a) droplet terminal velocity 
which depends on d (refs. 14 and 15) and on air properties, and (b) Sh, which 
depends on both d and V^, it is desirable to express the product Sh/(d^V^) 
in terms of d only to facilitate calculations of pH and A for the monodis- 
perse rainfall case. 

Values of Sh/(d2v^ were calculated for convective diffusion of HCS. (in 
air at 15^C and 0.85 atm) to water droplets falling at terminal velocity. 
Equation (4) was used to calculate Sh. The water droplet terminal velocities 
tabulated in Mason (ref. 15) for 20^C and 1 atm were corrected to conditions 
of 15^C and 0.85 atm by applying a factor of 1.0665. This factor was calculated 
from an expression developed by Foote and Du Toit (ref. 16), and it is assessed 
by Mason (ref. 15). Results which illustrate the slightly curvilinear power-law 
variation of Sh/(d^V^) with d are shown in figure 1. A best-fit straight 
line is described by 


= 0.022d-2*15 (13) 

d2v„ 

OO 


A more accurate fit, also shown in figure 1, is 


= 2.61 (10d)“2*15 + 0,175 (iQd) (14) 

d2v„ 

OO 


Either equation (13) or equation (14) can be substituted into equa- 
tions (7a), (9), and (12a) to obtain expressions that characterize washout by 
monodisperse rainfall. Use of the simpler but less accurate expression, equa- 
tion (13), results in 


pH = log 


1.125 X I07d2*15 
p(HCfi-)z 


p(HCfi,)z = Po(HCil)z exp (-6. 86 x 10~'^d"2. ISjjt) 
A = 6.86 X 10“‘^Hd“2.15 


(15) 


(16) 

(17) 
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where p(HCil)z is the effective vertical HCil(g) column density (in ppmv-m) 
t seconds after the onset of rain deposition. 

Substitution of equation (16) into equation (15) yields a combined expres- 
sion that can be used to calculate the pH of monodisperse- droplet rain at ground 
level any time after onset of rainfall deposition; 


1.125 X lO^d^.lS ^ o -.c 

pH = log + 2.98 x lo-'^-Z.lSHt (18) 

Po(HCil)z 


Note that applies only to H*** derived from HCt; also, t = 0 applies in a 
practical sense when the idealized monodisperse rainfall first reaches ground 
level. 


Washout Coefficient for Polydisperse Rain 

For naturally occurring polydisperse rain, droplet size spectra can be 
characterized by various frequency-distribution functions d(n)/d(d) (in 
drops/cm^), which depend on rainfall rate and type of rain (refs. 14, 15, 
and 17). These functions can be used to define (see ref. 14) an integrated 
intensity; by using equation (11), 


H 


. d(n) 

60007T d^V^ d(d) 


(19) 


and the washout coefficient for highly soluble gases, from equation (10), 


A = I d(Sh) d(d) (20) 

Jo ‘a(d) 


Note that the generalized integral expression for A in terms of rain-droplet 
collection efficiency E(a,d) is given (refs. 11 to 14) by 


A = 



d(n) 

E(a,d)V^ d(d) 

d(d) 


( 21 ) 
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where collection efficiency for convective diffusion with irreversible 
absorption is 


E(a,d) 


(Sh) 

CO 


( 22 ) 


Expressions for £(a,d), which apply for Brownian diffusion, interception, and 
impact collection of aerosol particles and/or droplets of radius a rain- 
drops of diameter d are discussed by Slinn (ref. 14) and Pruppacher and iG.ett 
(ref. 17). 

Although equation (20) can be integrated by using numerical substitution to 
obtain A at various H for any given raindrop spectra, the following alterna- 
tive approach was taken to eliminate the implicit dependence on (contained 

in Sh) . This approach greatly simplifies subsequent numerical integrations for 
A by having functions that contain only d and H. Values of d(Sh) were 
calculated for convective diffusion of HCA (in air at 15®C and 0.85 atm) to water 
droplets falling at terminal velocity. Results which illustrate the slightly 
curvilinear power-law variation of d(Sh) with d are shown in figure 2. A 
best-fit straight line is 

d(Sh) * 52dV3 (23) 


Now that d(Sh) is accurately expressed as a simple function of d only, 
substitution of equation (23) into equation (20) leads to a convenient integral 
expression which can be used to evaluate A for rains of known droplet size 
distribution 


A = 52 itD;^ 


r d5/3 
Jo d(d) 


d(d) 


(24) 


Washout Coefficient Based on Marshall-Palmer Distribution 

The Marshall-Palmer (M-P) raindrop size distribution (ref. 18) has been 
used in a number of studies to calculate precipitation scavenging characteris- 
tics. (See ref. 14.) The expression for the M-P differential number distribu- 
tion for raindrops with diameters between d and d + d(d) is given by 


d(n) 

= 0.08 e3q?(-41H“0*21d) 

d(d) 


14 


( 25 ) 



Before the M-P distribution is adopted for use in this analysis, it is 
important that its principal attributes and limitations be recognized. Mason's 
initial assessment (ref. 15) of the M-P expression for raindrop size distribu- 
tion, equation (25) , was that it provided a good average fit to the experimental 
data of Laws and Parsons (ref. 19), apart frcan a tendency to overestimate the 
numbers of small drops. Based on four more recent independent sets of data (up 
to 1962, ref. 20), Mason (ref. 15) concluded, "While the Marshall-Palmer 
distribution represents quite well the average size distribution of a large num- 
ber of samples from continuous cyclonic or warm-frontal rain, or the character- 
istics of the rain averaged over a considerable interval of time, samples taken 
over periods of only a few minutes may show considerable deviations frcMn these 
spectra." Slinn (ref. 14) discussed the problem of defining an appropriate 
raindrop size distribution, especially for predictions of the washout coeffi- 
cient. He compared empirical formulas for seven different size distributions 
and asserted, "The difference between real data and these smooth empirical func- 
tions can be very large, especially for convective storms where frequently 
they're not even qualitatively similar." Slinn also pointed out that even for 
continuous frontal rain, as studied by Waldvogel (ref. 21) and Joss, the param- 
eters in the M-P distribution had large variations during a single storm. 
Finally, Pruppacher and Klett (ref. 17) provided a very recent and comprehensive 
assessment of the M-P and other distribution functions used to characterize 
rainfall under a variety of conditions. 

By considering the limitations of the M-P raindrop size distribution and 
also recognizing that this distribution still remains an acceptable standard for 
predictive comparisons, it can now be applied to the determination of HCfi/(g) 
washout coefficient as a function of rainfall intensity. Substitution of equa- 
tion (25) into equation (24) with = 0.187 cmVsec followed by numerical 

integration over the finite droplet size range 0.01 S d ^ 0.6 cm for each of 
several H-values in the range 0.1 i H ^ 100 mm/hr leads to the results shown 
in figure 3. The calculated results (circles) are well characterized by the 
empirical fit 


A = 1.80 X IO^^hO-SSS 


(26) 


Washout Coefficient Based on Kelkar Raindrop Intensity Data 

Experimental data in the form of raindrop intensity (in drops/m^-sec) ver- 
sus size can also be used to deduce the washout coefficient and its dependence 
on rainfall rate. Kelkar (ref. 22) published the results of an extensive set 
of ground-level raindrop measurements and gave the intensity of raindrops as 
a function of droplet diameter for several different rainfall rates. The data 
used in the present analysis represent averages of grouped results (table 2 of 
ref. 22) obtained during 8 rain periods, with 49 000 drops counted. They are 
confined to general rains, as opposed to Kelkar 's other studies of drizzles, 
light showers, moderate showers, and monsoon thunderstorm rains (ref. 23) . The 
data are grouped into sets, with averaged values of the rainfall intensity and 
the number of drops/m^-sec tabulated for 21 intensity groups (0.20 to 34 mm/hr) 
and 16 diameter groups (0 to 0.25 mm up to 3.75 to 4.00 mm in 0.25-mm 
increments) . 


15 



In order to utilize the Kelkar raindrop intensity data to calculate A 
for polydisperse rains, we first denote these incremental data by An^ for the 
ith class of mean diameter di and corresponding terminal velocity V^, and 
then express the incremental data in terms of An^; 


Ani = lO^v^Ani 

Equation (20) is rewritten as 


(27) 


A = 



d(d) 


(28) 


Dividing equation (28) by equation (19) gives 



(29) 


When equation (13) is substituted as the simple power-law expression for Sh/d^v^,, 
the respective integrals of equation (29) are approximated by finite sums, 
and 10“^ Ani substituted for Ani, following expressions for A are 

obtained: 


N 


A = 


y (o.022di-2.15jdi3 Ani 

HDab 6l ^ 


6000 


N 

y di3 Ani 
i=l 


(30) 


Sh 


( 31 ) 
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where the weighted mean is defined by 



0.022 


N 

i=l 


85 


An^ 



( 32 ) 


Kelkar's averaged raindrop intensity data were used to calculate ^Sh/d^V^^ 
through application of equation (32) to his 21 sets of cOTibined droplet 'inten- 
sity data, each corresponding to an average rainfall rate. The results shown 
in figure 4 indicate that 


m 



3.47H"0-375 


(33) 


represents a good empirical fit of the data over the entire experimental range of 
rainfall rate. Substitution of equation (33) into equation (31) results in an 
expression for the washout coefficient for Hdl (g) based on D ;^3 = 0.187 cm^/sec 
and the Kelkar raindrop intensity data: 


A = 1.08 X 10“4h5/8 


(34) 


Respective values for A calculated directly from equation (30) for each of 
Kelkar *s 21 sets of averaged raindrop intensity data are shown in figure 3 
(squares) and are faired with a straight line representing equation (34). 

The somewhat more accurate expression for Sh/d^V^, equation (14), could 
have been substituted into equation (29) to define a refined (and more complex) 
finite sum expression for A analogous to equation (30) . But this treatment 
of the Kelkar data did not appear to be justified in view of the relatively 
large drop size interval used; e.g., especially over the inclusive lower bound 
interval (0 to 0.025 cm), where the relative difference in accuracy between 
equations (14) and (13) is the greatest. 

The expression for Hdl (g) washout coefficient based on Kelkar *s data, equa- 
tion (34), can be compared with Engelmann's graphically illustrated A for 
iodine (I 2 ) washout (ref. 11), which was based on an independent numerical eval- 
uation of equations (3) and (20) using the same Kelkar rain intensity data. The 
present author's measurement of Engelmann's slope (0.62) for log A versus 
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log H indicated nearly exact agreement with the present H^/S dependence. 
Also, the estimated ratio of multiplying factors (2.6 for 25*^ and 1 atm) was 
very close to an estimated ratio of molecular diffusion coefficients, as 
calculated frcm the inverse square root of the molecular-weight ratio (2.64). 
Thus, the present determination of A for HCil(g) washout, using Kelkar's data, 
is consistent with Engelmann's independent prediction of A for washout of 
I 2 vapor. 


HC?-(g + aq) Washout From SRM Exhaust Scavenging 

The HC<l(g) washout coefficients can now be compared with scavenging results 
obtained from experimental studies of solid-propellant exhaust. Consecutive 
studies by Knutson and Fenton (ref. 24) and Fenton and Purcell (ref. 25) of the 
Illinois Institute of Technology Research Institute (IITRI) are of particular 
relevance, since they represent the only published attempts known to this writer 
to measure rain scavenging rates in laboratory studies of solid rocket exhaust. 
Because of the potential importance of results derived fran these data, these 
studies are first thoroughly reviewed and critiqued in appendix A. Then 
the Fenton-Purcell (F-P) result is corrected, and the data are reanalyzed in the 
following two subsections. 


Correction of Fenton-Purcell A Expression 

Close inspection of the A expressions derived by Fenton and Purcell 
(see appendix A and also p. 53 of ref. 25) clearly indicates that both are 
incorrect algebraically. The ratio written to be 

instead of This error has a relatively large effect on the calcu- 

lated values of A at all relevant HCH concentrations (0.1 to 100 ppmv) . The 
power-corrected F-P equation should be 

A = 4.21 X 10-8pgirMA"°*^^^H0*773 (35a) 


when Pair expressed in g air/m^. 

Substitution for Pair ^A power-corrected F-P A expression 

leads to 


A = 1.41 X 10"4p(HCd) ■0*''76h0.773 


(35b) 


for the presently selected atmospheric conditions at 15°C and 0.85 atm. Results 
calculated from this power-corrected expression are shown in figure 5 for HCH 
concentrations ranging frcsn 0.1 to 100 ppmv. Comparable results for the 
original F-P expression (not shown) were incorrectly low by a factor of 10 for 
1 ppmv HCSl, for example. Finally, note that the single solid line representing 
this writer's analysis of the F-P data (shown later) crosses through the various 
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power- corrected F-P results for different HCfi, concentrations. Details of the 
analysis are derived in the next subsection of this paper. 

Although the HCH concentration dependence of the power- corrected F-P ej^res- 
sion is relatively weak, it is evident in figure 5 that significant increases in 
A with decreasing p(HCil) would result for a realistic range of p(HCil). A 
review of the F-P results, however, indicates that the statistical significance 
of the derived HCil concentration dependence was not demonstrated. Furthermore, 
this dependence appears doubtful because of the large data scatter, inclusion of 
the previously discussed high- concentration reference 24 data, and the presence 
of large run- to-run variations in relative humidity. 


Alternate Analysis of Fenton-Purcell Scavenging Data 

An alternate analysis of the F-P data is now presented to (a) examine the 
validity of the corrected washout coefficients and (b) find a simpler expression 
for A which is consistent with the data. Note that Fenton and Purcell based 
their reference 25 analysis on the major a priori assumption that the parameter 
^A/^^air^oo^^ correctly scales the combined effects of variable droplet diameter 
and droplet fall speed on scavenging rate. In order to (a) allow a straightfor- 
ward examination of this assumption and (b) eliminate the apparently weak depen- 
dence of A on HCil concentration (statistical significance uncertain), it was 
assumed provisionally that the experimental rate of chloride scavenging is 

directly proportional to chamber HCil (g + aq) concentration C(HCil) for each of 
the three droplet-size — fall-speed combinations investigated. 

The present author's calculated (from tabulated data) scavenging rates 
are plotted against C(HC^l) in figures 6, 7, and 8 for average droplet diam- 
eters of 0.055, 0.11, and 0.30 cm, respectively. All 29 usable data runs, con- 
tained in table VI of reference 25 and asserted to be reliable, were included in 

these log-log plots. Experimental chamber relative humidities are shown adja- 
cent to each point. Although the data scatter is substantial and some runs are 
clearly of lower quality than others because of the relatively large rainwater 
chloride corrections required sometimes (starred data), the entire data set is 
used throughout the present analysis. The 45^ lines in each plot represent mean 
values of R;^/C(HC^) in accordance with the previous assumption of linear con- 
centration dependence. Inspection of the data fails to indicate the presence of 

a statistically significant effect of chamber relative humidity on which is 

in agreement with the F-P assessment. 

The mean values of R^/C(HC?.), which effectively represent washout coeffi- 
cients on a per-droplet basis for each monodisperse droplet size, are plotted as 
a function of d in figure 9. The least-squares expression 

Ra/C(HC«.) * 1.24 X I0“7d2.13 

is a good correlation of the data. The vertical bars indicate 1 standard devia- 
tion of the mean values. Since the molar collection rate per droplet is given 
by = R;^/36.45 (prime denotes HCS- (g + aq) collection), equation (36) is 
equivalent to 
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( 37 ) 


10®Wa 

= 94.5d2-13 

CfC (HCH) 


when the present value for Cf at 15*^ and 0.85 atm is used. 

The corresponding washout coefficient for polydisperse rain is now defined 
by an integral expression, similar to that defined previously, as 


A = 


/ 


d(n) 

CfC(HCfi.) d(d) 


d(d) 


(38) 


After substituting the M-P raindrop size distribution function, equation (25) , 
A can be evaluated numerically from 


A = 


r 


6 

01 


(94.5d2.13) [0.08 ejg) (-41H“0* 21d) ] d (d) 


(39) 


Calculated values for A at various rainfall rates (circles) are shown in 
figure 10. An excellent fit of these results for HCA (g + aq) washout is 


A = 1,52 X 10“'*h0*658 


(40) 


Equation (26) for UCl {g) washout is also plotted in figure 10. Note that 
both apply to the same ambient and droplet size distribution and integration 
conditions. Surprisingly, the respective values of A for these independent 
cases are remarkably close over the normal range of rainfall rate. In fact, 
they are statistically identical when the data scatter in figure 9 (±30 per- 
cent) is considered. Unfortunately, the reason for this coincidence of 
scavenging results cannot be evaluated quantitatively without measurements 
of aerosol size distribution and analytic considerations of aqueous aerosol 
growth and scavenging, as discussed previously. 


Summary of A Expressions for HCil(g) Washout 

The HCl washout coefficients based on the M-P raindrop size distribution 
and the Kelkar raindrop intensity data are characterized by a generalized 
empirical expression 

A = AH*^ (41) 
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The values of A and b for these HCt(g) washout cases are summarized in 
table II, along with values for the geometric mean of the two determina- 
tions. The author's previously evaluated expression (ref. 3), based on an 
alternate treatment of the same Kelkar data, an outdated binary diffusion coef- 
ficient (0.20 cmVsec), and standard air properties at 25^C and 1 atm, is shown 
for comparison; although the previous value of A, 1.11 x 10”^, was very close 
to the present one, 1.08 x 10“^, this small difference partly reflects compen- 
sating effects. 


TABI£ II.- SUMMARY OF BCl WASHOUT OOEFFICIENTS IN TERMS OF GENERALIZED EXPRESSION A == AH^ 


HCJU ) 

Temperature, 

OC 

Pressure, 

atm 

Source 

Scavenging 

rate 

Raindrop 
distribution; 
A- equation 

A, 

sec“^ 

b 

g 

15 

0.85 

Present 

Modified 

Frossling 

M-P; eq. (26) 

1.80 X 10-4 

0.565 

9 

15 

0.85 

Present 

Modified 

Frossling 

Kelkar; 
eq. (34) 

1.08 X 10^4 

0.625 

g 

15 

0.85 

Present 

Modif ied 
Frossling 

Geometric 
mean of M-P, 
Kel kar ; 
eq. (42) 

1.39 X 10-4 

0.595 

g 

25 

1.00 

Previous 
(ref. 3) 

Modified 

Frossling 

Kelkar 

1.11 X 10“4 

0.625 

g + aq 

15 

1.00 

Knutson 

and 

Fenton 
(ref. 24) 

Data with 
50% of 
modified 
Frossling 

M-P 

.83 X 10“4 

0.567 

g + aq 

15 

0.85 

Presently 

corrected 

F-P 

r es ul t , 
eq. (35a) 

Data 

M-P; eq. (35) 
(note ppmv 
Ed in 
A- factor) 

1.36 X 10”5 p(HC«.)-0.176 

0.773 

g + aq 

15 

0.85 

Present 
analysis 
of F-P 
(ref. 25) 

Calculated 

from 

data 

M-P; eq. (40) 

1.52 X 10“4 

0.658 

I 


The SRM exhaust rain scavenging results reviewed on pages 18 to 20 and 
41 to 44 of this paper are also summarized in table II. While Knutson and 
Fenton (ref. 24) used 0.173 cm^/sec for their HC0/(g) diffusion coefficient at 
15^C and 1 atm and air properties under the same conditions, multiplication 
of their A value 2 (= 1.66 x 10“^) yields a value of A which is approxi- 
mately equivalent to the present result (first row) when the difference in 
atmospheric pressure is fully accounted for. Although the corrected F-P 
expression (eq. (35)) yields A values that agree roughly with those derived 
from the reanalysis of their data (eq. (40)), the latter is preferred because 
of its simplicity and also the absence of a clear e3q>erimental dependence of 
A on HCi!/ concentration and humidity. 
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Two recommendations are now offered. First, since (a) equation (40) for 
HCJl (g + aq) washout is statistically equivalent to equation (26) for HCA (g) 
washout over the relevant H range and (b) the physics of HCSl(g) washout are 
more firmly established, it is provisionally recommended that rocket exhaust 
HC<1 washout analyses be based on values of A which apply to rain scavenging 
of HCS< (g) . If at scnie later point new aerosol size distribution information 
and detailed SRM cloud micro£*»ysics and scavenging calculations lead to well- 
defined alternative expressions that apply, for example, at high humidities 
(approximately ^ 95 percent), then reconsideration will be necessary. Second, 
it is recommended that the geometric mean of the M-P and Kelkar integrated 
HCfi. (g) washout coefficients (both based on the modified Frossling equation) be 
used for predictive purposes; i.e.. 


A = 1.39 X io-4hO*595 


(42) 


This compromise is considered both a reasonable and conservative expression for 
characterizing HCH (g ) /HCS/ (g + aq) washout under the presently selected atmo- 
spheric conditions at 15°C and 0.85 atm, since neither of the raindrop size 
intensity distributions is considered unequivocal. 


Resultant pH of Polydisperse Rain 


Using the HCS,(g) washout coefficient, 
ring, polydisperse rains by equation (41), 
pH of rain at ground level can be derived, 
case, equation (12a) is rearranged as 


defined for average, naturally occur- 
A = AH*^, a general expression for the 
First, for the monodisperse rain 


Sh _ 6000A 
d2v„ ®AB 

OO 


and then substituted for Sh/(d2v^) in equation (7a) to eliminate the explicit 
dependence of pH on d, and Sh; thus 


pH = log 


H 

3600CfAp(HCA)z 


(43) 


Substituting equation (9) for p(HC8,)z gives 


pH = log 


H 

3600CfApo(HC2,)z 


At 

+ 

In 10 


(44) 
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Note that for ground deposition considerations t > 0 pertains to time lapse 
after uniform-size raindrops of diameter d (used in eg. (18)) first reach 
ground level. Thus, while washout actually occurs at seme earlier average time 
t - t', its effect is not measured at ground level until t > 0. 

Since A is the only term in equations (43) and (44) which depends on the 
distribution of raindrop size (assuming fixed H), equation (41) may be substi- 
tuted for A to yield an ejq)ression for the £« of polydisperse rain at ground 
level; thus 


pH = log 


Hl-b 

3600cfApo(HCA)z 


AHt»t 

+ 

In 10 


(45) 


When the geometric mean A in equation (42) (i.e,, A = 1.39 x 1 0“^ and 

b = 0.595) is used, equation (45) reduces to 


pH = log 


5.55 X IO^hO- 405 
Po (HC5.)z 


+ 6.04 X 


(46) 


Several features of equation (45) should be noted. First, it consists of 
two A- and H-dependent terms: one relates to HCt(g) column density at the 

onset of rainfall deposition, and one represents the effect of progressive wash- 
out. Second, it represents the mixing-cup (volume-weighted average) p*[ of steady 
rainfall received at ground level; thus, it has been assumed that differences in 
transit times 6t* for the various droplet sizes which contribute most to 
(dominate) volumetric average pB are small compared with incremental values of 
t which alter pH significantly in equation (45) . In practice, an instantaneous 
mixing-cup average pH would not be realized exactly, since droplets of different 
size do have different transit times through the SRM cloud- Thus, a finite 
sampling time is appropriate for practical application of equation (45) , partic- 
ularly during initial stages of rainfall and, in general, for any natural rain- 
fall, since large spatial and temporal variations in raindrop size distribution 
usually occur during a given rainfall event. Note that for calculation of the 
final total acidity pHf of a rain that has an initial acidity pHq, the relation 


pHf = -log (10"£® + 10 


(47) 


can be used. This satisfies the problem of an inoorrect asymptotic limit for 
pH in equations (7), (15), (18), and (43) to (46) for very long washout times 
and/or very low vertical HCll(g) column densities. 
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SRM EXHAUST CLOUD DISPERSION 


Specific applications of the NASA/MSFC cloud rise preprocessor and the var- 
iously modified Gaussian multilayer diffusion models (MDM) of exhaust cloud dis- 
persion for tropospheric air quality predictions at ground level are described 
in references 26 to 32. Sane of these references focus on the inherent mathe- 
matical and physical limitations of this basically statistical analytic approach 
(ref. 32)/ parametric studies of the principal variables upon which the MDM 
models are based (refs. 26 and 32) / and comparisons of calculated rocket effluent 
concentrations at ground level with those obtained from other advective diffu- 
sion models (ref. 32) . 

The parameters used in this paper to characterize HCA source strength mo 
and vertical HC?/ column density a are documented by this author in refer- 
ence 6 and an unpublished paper. The original SRM cloud dispersion calculations, 
fron which these parameters were derived, were obtained through use of the NASA/ 
MSEC cloud rise preprocessor and Gaussian MDM described and documented in refer- 
ences 28 to 30 for operational prediction of toxic fuel hazards. Dispersion 
calculations using Model 4, version II of the MDM codes, designated here as 
MDM-4(II), comprise the basis of parameters used throughout this study. Simul- 
taneous inclusion of the precipitation scavenging subroutine is indicated by the 
designation of MDM-5(II). 


Standard Meteorological Regimes for Cape Canaveral 

The meteorological regimes selected for this work consist of the set of 
standard meteorologies (ref. 5) that were originally used to assess atmospheric 
dispersion of rocket exhaust effluents in the Cape Canaveral, Florida, area. 

They represent the major meteorological conditions which lead to overland trans- 
port and are likely to be encountered in the Cape Canaveral area. The meteoro- 
logical profiles of temperature, wind speed, and wind direction for all seven 
meteorological regimes are shown in figures 11(a) to 11(g), as obtained from 
references 33 and 28 (for sea breeze (SB)), and are discussed in appendix B. 
Meteorological data for specific launch occasions are generally tabulated in 
separate publications (e.g., refs. 34 to 36). 


Characteristics of Altitude-Stabilized SRM Source Clouds 

The SRM cloud rise preprocessor portion of the MDM-4(II) model numerically 
calculates the vertical and horizontal distribution of HCt which results from 
buoyant cloud rise with turbulent air entrainment under a specified potential 
temperature gradient (refs. 26 to 36). The resulting vertical source-strength 
distribution is then assigned a multilayer structure, based on the meteorologi- 
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cal profiles^ which includes a specific geometric shape in the vertical cross 
section to define the initial distribution of exhaust mass in the altitude- 
stabilized SRM cloud. In the case of Model 4, a conical body of revolution, 
having a diamond-shaped cross section symmetrical about the vertical centroid 
axis, is used to bound the horizontal layers. Examples for the SB, cold front 
passage (CFP) , and post-cold front (Post-CF) cases are shown in figure 12 and 
are discussed in appendix C. 


Source-Strength Considerations 

Although the vertical-line — source-strength input function (propellant mass 
burn rate versus altitude) is accurately known for any prescribed Titan III or 
Shuttle mission, the cloud rise process determines resultant vertical and hori- 
zontal SRM exhaust distributions, which apply after altitude stabilization 
has occurred. Thus, the calculated layer source strengths, which lead to corre- 
sponding initial HCA concentrations in the MDM subdivided stabilized SRM cloud, 
always differed from the accurately known input source-strength increments. 
Moreover, the resultant cumulative HCi source strength in the surface mixing 
layer was always smaller than the cumulative input of HCil mass exhausted up to 
the defined stabilization altitude. Finally, both the resultant layer and cumu- 
lative HCjl source strengths varied significantly as functions of the various 
potential temperature profiles. 

Input histories and resultant HCJl source strengths which pertain to the 
chemically similar SRM boosters for both Space Shuttle and Titan III launch 
vehicles are shown in figure 13. The cumulative HCA inputs to the MDM for both 
Shuttle (ref. 37) and Titan III (ref. 28) launches are based on propellant con- 
sumption rates and launch mission trajectories. They differ by a factor of 2.44 
up to an altitude of 2.0 km. The resultant preprocessed fall fair weather (FFW) 
cumulative HC?. source strength in the MDM stabilized cloud for Titan III is also 
shown in figure 13 for altitudes up to 2 km. The so-called original spring 
fair weather (SFW) cumulative source strength in the MDM stabilized cloud, used 
in previous acid rain calculations (ref. 3) , is shown as a single point (4 km, 

89 t HCA) . In view of the Mission 2, November 1973 design (ref. 37), it is too 
large by about 30 t for a 4-km upper cloud altitude boundary. 


Summary of Exhaust Cloud Characteristics 

Table III summarizes, for Titan III launches and the seven standard meteo- 
rological regimes at Cape Canaveral, Florida, SRM exhaust cloud characteristics 
deduced from application of the MDM-4(II) in reference 6 and an unpublished 
paper by G. L. Pellett. In some cases the mean transport wind speed for the 
cloud Uq, which relates distance to time in the model, differs significantly 
from the wind speeds of the individual layers. 
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TABLE III.- SRM EXHAUST CLOUD CHARACTERISTICS, AS DEFINED FROM APPLICATION OF M)M-4(II) TO TITAN III 
LAUNCHES, FOR SEVEN STANDARD METEOROLOGICAL CASES AT CAPE CANAVERAL, FLORIDA 


Meteorological regime 

Meteorology 

abbreviation 

Cloud - 
centroid 
height, 

Zjn, m 

Surface 
mixing- layer 
depth, 
m 

MDM upper 
Cloud 
boundary, 
m 

Effective 
cloud thickness 
at X ^ 5 tan, 
m 

Mean wind 
speed for 
cloud, 
Uq, id/s 

Mean wind 
direction 
(North = 0°), 
deg 

Cloud 

stabilization 

time, 

sec 

Calculated 
distcmce, 
pad to cloud 
stabilization 
(from Uc) , 
km 

Empirical 
virtual source 
distance, 

Xq, km 

HCH 
source 
strength, 
®o» 9 

Fall fair weather 

FEW 

1322 

1000 

2000 

1000 

6.127 

105 

364 

2.23 

10 

14.88 X io6 

Spring fair weather 

SFW 

1380 

2000 

2000 

2000 

7.0314 

140 

416 

2.93 

4.0 

14.48 

Low-level sea breeze 

LLSB 

959 

300 

2000 

1700 

6.752 

145 

199 

1.34 

3.5 

8.634 

Sea breeze 

SB 

993 

800 

2200 

1400 

9.923 

165 

247 

2.45 

5.0 

10.12 

Fair weather, pre-oold front 

FW, Pre-CF 

875 

200 

1400 

1180 

2.453 

224 

160 

.39 

3.0 

11.89 

Cold front passage 

CFP 

1230 

2000 

2000 

2000 

6.696 

66 

301 

2.02 

4.0 

15.07 

Post- cold front 

Post-CF 

1341 

1400 

2000 

600 

8.712 

69 

442 

3.85 

2.5 

11.30 

Post-cold front (pad abort) 

Post-CF 

1506 

1400 

2000 

600 

9.284 

69 

342 

3.18 

2.5 

15.86 


(pad abort) 
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The cloud stabilization times from the model were used to estimate downwind 
drift distances frcm the launch pad to the point where cloud altitude stabi- 
lization theoretically occurs. These cloud stabilization drift distances cure 
compared with the empirically determined virtual source distances Xq (Xq was 
needed to obtain straight-line power- law fits of a versus X). Finally, the 
resultant HCil source strengths for the stabilized SRM clouds are also sum- 

marized in table III. They apply from the Earth's surface to the MDM upper 
cloud boundary and are conserved throughout the SRM cloud dispersion history (in 
the absence of rain), since no loss terms were applied at the SRM cloud 
boundaries. 


Vertical HCt Column Density Determinations 

The MDM-4(II) calculated downwind vertical HCS. profiles and the ompirical 
power-law fits of a versus downwind distance X = X^g + Xq are illustrated in 
reference 6 and an unpublished paper by G. L. Pellett for each of the Titan III 
cases studied. The original SFW Shuttle case used in reference 3 is reproduced 
in figure 14, and the present power-law-decay-modified SFW Shuttle cases (two) 
are shown in figure 15. The solid line in figure 15 represents the original 
source strength; the dashed line, the present fully modified source-strength 
case. 


The decay expressions for vertical HCd column density a = aX“^ are sum- 
marized in table IV for the Titan III and Shuttle cases, along with HCil source 
strength, mean transport wind speed, and other relevant exhaust cloud properties 
derived frcxn the MDM-4(II) calculations. Corresponding time- dependent decay 
expressions for a are also given. The empirical fits for a were excellent 
throughout the data range, 1 ^ X i 100 km, for all the Titan III meteorological 
cases studied. It will be shown later that the demonstration of adequate power- 
law ejqjressions for a leads to substantial analytic simplification, especially 
for subsequent calculation and parameterization of acid rain characteristics. 
Some effects of meteorology on expressions for a in table IV are now dis- 
cussed. Although the deduced values of a can be referred to as an HC?. column 
density source strength at unit distance, the reader should recognize that a 
is inherently sensitive to Xq. In turn, Xq is influenced in the MDM-4(II) 
cloud rise calculations by the vertical profiles of temperature and wind speed 
as well as the placenent of vertical layer boundaries. Thus, the physical sig- 
nificance of a is approximate and is not precisely defined. 

A noteworthy internal consistency exists for the two Post-CFP meteorologi- 
cal cases, which apply to a normal launch and a pad abort situation. The 
respective Xq and 3 values were identical, and the ratio of a values (and 
a values) (0.71) was identical to the ratio of values. In the absence of 

other similar paired calculations for Titan Ill's, it is not known whether a 
simple linear scaling law would apply to other meteorological cases. 

A roughly similar comparison can be made between the fully modified Shuttle 
SFW and Titan III SFW cases, but it is not as exact since additional multilayer 
structure bounded at an altitude of 4000 m was used for calculating the Shuttle 
case, whereas the Titan III case was terminated at 2000 m. While the respective 
3 values differed (1.64 and 1.98), the Xq values were the same (4.0 km), and 
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TABLE IV.- EMPIRICAL SRM EXHAUST CLOUD PROPERTIES CHARACTERIZING DOVWWIND DISPERSIVE DECAY OF VERTICAL HCH OOLUW DENSITY 

Deduced frcm application of MDM-4(II) to Titan III launches for seven standard meteorological cases at Cape Canaveral, 
Florida; the present power- law- decay- modified SFW Shuttle case is shown for both originally used source strength 
ihq * 89 X 10® g HCi (ref. 3) and presently reduced source strength m^ * 61 x 10® g HCil consistent with more recent 
Shuttle mission design (ref. 37) 


Meteorology abbreviation 

Cloud- 
centroid 
height, 
Znif m 



M)M 

upper cloud 
boundary, 
m 

Effective 

cloud 

thickness at 
Xcs ^ 5 tan, 
m 

Mean wind 
speed for 
cloud, 
Uq, m/8 

HCil 
source 
strength, 
mo, 10® g 

Empirical peurameters for vertical HCi col min density, a, ppinv-m 

G * 0t(1^CS 


o«-B 

CJ ~ 01(^05 

to)-B - 

ext- 6 

Ot, 

ppmv^m 

Xq. 

km 

a 

Of 

ppmv-m 

to» 

hr 

3 

1 

Titan III cases 

FFW 

1322 

2000 

1000 

6.127 

14.88 

1.6 X 10® 

10 

0.84 

11.9 X 1q3 

0.453 

0.84 

SFW 

1380 

2000 

2000 

7.0314 

14.48 

2,6 

4.0 

1.98 

.433 

.158 

1.98 

LI£B 

959 

2000 

1700 

6.752 

8.634 

4.3 

3.5 

1.32 

6.35 

.144 

1.32 

SB 

993 

2200 

1400 

9.923 

10.12 

14.0 

5.0 

1.81 

2.17 

.140 

1,81 

PW, Pre-CF 

875 

1400 

1180 

2.453 

11.89 

2.0 

3.0 

.93 

26.3 

.340 

.93 

CFP 

1230 

2000 

2000 

6.696 

15.07 

4.5 

4.0 

1.93 

.968 

.166 

1.93 

Post-CF 

1341 

2000 

600 

0.712 

11.30 

.60 

2.5 

1.13 

1.22 

.080 

1.13 

Post-CF (pad abort) 

1506 

2000 

600 

; 9.284 

1 15.86 

.85 

2.5 

1.13 

1.65 

.075 

1.13 


Space Shuttle cases 


SFW 

(original source 

strength) 

— 

4000 



7.50 

89.0 

— 

14.6 

4.0 

— 

1.64 

6.56 

0.148 

1.64 

|SFW 

(modified source 

strength) 

— 

4000 

— 

7.50 

61.0 

10.0 

4.0 

1.64 

4.49 

.148 

1.64 
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the ratio of a values (3.85) was fairly close to the ratio of niQ values 
(4.21). In terms of corrected distance from launch site, the point at which the 
ratios of a values and m^ values were equal to 4.21 was X = 1.3 km, which 
is still close to the launch site. Thus, these comparisons between Titan III 
and Shuttle for SFW meteorology suggest that approximately linear scaling of CJ 
with source strength applied. 

WASHOUT OP DISPERSING SRM EXHAUST CLOUDS 
Application to a Dispersing SRM Cloud 

Equation (45) and related e3q>ressions for predicting rain pH. and ground 
deposition of can be applied to independently dispersing SRM clouds once the 
temporal and spatial variation of vertical HCH column density is specified in 
terms of time and/or ground coordinates. The power-law decay of o which 
resulted frcxn application of MDM-4(II) is now defined by the integral expression 


Po^ 2 (HC! 1 ) dz = Po(Ha)z = aX-S 


where Pq HCZ {g) concentration at altitude z, Pq(HC^) repre- 

sents a vertically averaged concentration over z, X is the downwind distance 
frcxn launch site, and a and 6 represent onpirical constants. 

For simplicity, it is now assumed that the SRM cloud geometry is an expand- 
ing right circular cylinder of diameter Dx (axis vertical) with mass conserva- 
tion for HC?, at concentration Pq(HC< 1). Since the scavenging model applies 
to irreversible absorption and thus exhibits a direct proportionality between 
rain molarity and HC<l(g) column density, the actual distribution of HC^(g) in 
the vertical direction is unimportant. In order to calculate rain and ground 
deposition per unit area along the SRM cloud-centroid path (CC path), the hori- 
zontal distribution of vertical HC?/ column density along the projected cloud 
diameter is normally assumed uniform. For calculation of two-dimensional acid 
footprints, or ground deposition isopleths, the variation of vertical HC^ column 
density with ground coordinates X and Y must either be uniform beneath the 
SRM cloud or else specified. 

Substitution of equations (41) and (48) into equation (9) yields an expres- 
sion for the reduction of vertical HC?/ column density by dispersive decay and 
^ washout for X > 0 and t ^ 0; thus 

p(HCH)z = ax"^ e^(-AHt»t) (49a) 
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Substitution of equation (48) is made with the tacit assumption that SRM cloud 
dispersion is unaffected by the development and presence of an overriding rain, 
thus, the SRM cloud is denoted as independently dispersing. 


Now let Uq (in m/sec) represent the steady mean wind speed that character- 
izes SRM cloud transport, and Xj. (in km) represent the closest distance from 
launch site where a steady, overriding, vertical rainfall commences before 
arrival of the SRM cloud. Thus t = 1000(X - Xr)/!!^ sec represents the effec- 
tive time lapse after the SRM exhaust cloud encounters rainfall. Substitution 
of this expression for t in equation (49a) leads to an equation, in terms of 
downwind distances only, which describes the combined effects of dispersion and 
washout on vertical HCA column density for 0 < Xj. ^ X; thus 


p (HC<1 ) z = ax~^ 


exp 


-1000Ah‘^(X - Xr)/Uc 


(49b) 


Resultant pH From a Dispersing SRM Cloud 

An expression for the pH of polydisperse rain at ground level is obtained by 
substituting equation (48) into equation (45) and letting t = 1000 (X - Xj.)/Uq. 
For X ^ Xj. ^ 1 km. 


pH = log 


H 


1-b 


3600cf Aa i 


+ B log (X) + 


AH^ 


/X - Xr 


In 10/\Uc/1000i 


(50) 


For the geometric mean A (eq. (42)), equation (50) reduces to 


pH = log 


5.55 X 10^h0*405 


a 


+ B log (X) + 6.04 X to-2hO*595 


X - X, 


u. 


(51) 


In equation (50) there are three terms which affect rainwater pH; The 
first relates to a combination of rainfall intensity, washout coefficient, and 
HCS- source strength at unit distance; the second depends solely on SRM cloud 
dispersion; and the third accounts for progressive removal of HCA (g) through 
washout. Exclusion of the latter term, which is equivalent to setting X = X^, 
allows formation of an expression for potential rain pH, defined herein as the 
rainwater pH at ground level just after the first onset of rain. Succinctly, it 
represents the maximum rain acidity to be expected in the event of rain occur- 
rence downwind; thus 


PHpot = log 


H 


1-b 


3600cfA « I 


+ B log (Xr) 


(52) 
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Resultant Acid Deposition From a Dispersing SRM Cloud 


The areal deposition rate (flux) for aqueous at any point X on the 
SRM cloud-centroid (CC) path where washout is occurring overhead, is (averaged 
over passage of finite SRM cloud) 


Gx = 36.45 X 


(53) 


For a uniformly mixed cylindrical SRM cloud of effective diameter conserva- 

tion of HCS/ (in the absence of wet or dry deposition, absorption at the Earth's 
surface, and convective losses at the top) dictates that 


mo = 36.45Cf ^TTDxV4)(aX"^) 


(54) 


Thus the effective SRM cloud diameter along the CC path is 


Dx = 


m. 




36.457TCf a 


1/2 


(55) 


Since the time lapse for complete passage of the SRM cloud over any point on the 
CC path is D^/U^, the stationary ground receiver on this path will experience a 
total deposition per unit area of 


Gx = GxDx/3600Uc 


(56) 


Substitution of equations (53) and (55) into equation (56) yields an intermedi- 
ate result in terms of rainwater pH; thus 


3600/\TTCfy \^^2j 


(57) 


Substitution of equation (50) for pH and rearrangement of terms leads to a gen- 
eralized predictive expression for cumulative (total) areal HCS/ deposition along 
the CC path for all X ^ Xj-: 


^x “ 


-(36.45Cf) 

7T 


1/2 


1 /2 

expjj-1 000 AH*^q”^ (X - Xj-)J (58) 
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For the geometric mean A (eq. 42)), equation (58) reduces to 


Gx = 5.68 X 10-6h0*595d^-1 ^nioaX-3y/2 exp J^-0. 1 39r 0- 595u^,-l (x - Xj-)j (59) 

Ground receivers located at some distance Y (in km) normal to the CC path will 
experience, for the well-mixed cylindrical-cloud case, the same transient pH 
value but cumulative areal depositions of 


’x,y 


= 2Gv 


0.25 


1/2 



(60) 


for all Y < Dx/2. 

The total HC<1 deposition experienced at X = Xj. after the first onset of 
rainfall and for complete SRM cloud passage is defined herein as potential 
deposition and is obtained by setting X = Xj. in equation (58). For the pres- 
ent value of Cf (3.60 x 10“^ mol air/an^), the potential deposition is simply 

Gx,pot = 0.0409(AHbUc"^)(moaX-3y/2 


Collective Results for Dispersive Decay of a and Potential Rain pH 

The power-law decay of o for eight Titan III cases is shown as a function 
of corrected distance from launch site X in figure 16 and also as a function 
of elapsed time after launch X/3.6 Uq in figure 17. The fully modified SFW 
Shuttle case is shown as a dashed line in each figure. Two additional ordinate 
scales, illustrating the corresponding values of potential rain pH for rainfall 
rates of 1.5 and 25 mm/hr, are also shown in figures 16 and 17. This dual scale 
is deduced from either equation (46) with t = 0 or equation (52) using the 
geometric mean washout coefficient. 

The most noticeable feature of figures 16 and 17 is that the dispersive 
decay of a and potential rain pH differ greatly among the seven standard 
meteorological regimes. A range of more than 2 orders of magnitude in a and 
potential acidity (more than 2 pH units) is spanned at X ^ 100 km downwind 
and/or t ^ 2 hr postlaunch elapsed time. At shorter distances and times the 
total span in a and potential acidity is somewhat less but still exceeds an 
order of magnitude for X > 10 km and t > 0.2 hr. These results tend to con- 
firm earlier expectations of large variations in atmospheric dispersion under 
widely different meteorological conditions, and they serve to emphasize the need 
for developing a comprehensive atmospheric dynamics model (e.g., see refs. 10 
and 38) to deal effectively with this large source of variability. 

A second important feature of figures 16 and 17 is that environmentally 
significant potential rain acidities (e.g., pH ^ 1.5 results in short-term 
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damage to vegetation; see refs. 39 and 40) are shown to be possible for rela- 
tively large downwind distances from launch site and for long elapsed times. 

If the validity of the two meteorological cases, which tend to define upper 
bounds for a (ffw and FW, Pre-CF) , is accepted, then initial rain pH values 
between 1.0 and 1.5 (for 1.5 ^ H ^ 25 mm/hr) could occur at downwind dis- 
tances up to 50 km and elapsed times up to 5 hours. Although these estimates 
are derived from the two least dispersive meteorological cases, they are not 
considered excessively conservative. For example, somewhat more severe stagna- 
tion conditions are conceivable; no additional uncertainty bounds have been 
imposed; and the real possibility^ of a significantly increased raindrop collec- 
tion efficiency due to a predominance of large {^10-um diameter) aerosol drop- 
lets has not been accounted for. While the predictions in figures 16 and 17 
also indicate values of potential rain pH of 2 or less at downwind distances up 
to 200 km and elapsed times exceeding 10 hr, one must recognize that the valid- 
ity of these longer range results becomes less certain with increasing distance 
and time. 

Scxne competing effects, which both reduce and increase the potential rain 
acidity, are now identified. First, the MDM calculations do not include loss 
of HCt from the SRM cloud's upper boundary or sorption at ground level. 

Also, they do not include variable advection (wind speed and direction) for the 
respective MDM cloud layers. Inclusion of these processes would tend to reduce 
a and, hence, increase potential rain pH at progressively larger distances. 
Second, the variances used in the present MDM-4(II) calculations were origi- 
nally based on field data which applied to relatively small-scale atmospheric 
turbulence (refs. 34 to 36). It is well known (ref. 42), however, that as the 
relevant turbulence scale increases for atmospheric dispersion of very large 
plumes, the appropriate variances and overall dispersion rates tend to become 
weaker power-law functions of the characteristic scale size (e.g., source cloud 
diameter) and downwind distance. Thus, the exponent 3 , which characterizes the 
decay of a with X, could decrease to «0,5 at large X, in a fashion similar 
to that observed in large-scale dispersion studies (ref. 42) and that actually 
observed for Titan III exhaust clouds. The most extreme and best documented 
example out of eight Titan III cases observed thus far (ref- 43) was for the 
September 5, 1977, launch. Here, 45 data points for maximum in-cloud HC^ con- 
centration, obtained over the postlaunch period of 3 to 300 min, were character- 
ized by 3 = 0.54 in the empirical expression 


C(HCil) = 55t~0*54 


In summary, the sets of potential rain acidity results in figures 16 and 17 
characterize the range of dispersive decay O and initial rain (potential) pH 
that might be expected for Titan III launches. Respective envelopes bounding 
the collective results would encompass all seven standard meteorological condi- 
tions at Cape Canaveral but would not include the other uncertainties mentioned. 


^Informal communication of unpublished results obtained by A. N. Dingle 
and B. Heikes, University of Michigan, Ann Arbor, Michigan. See also refer- 
ences 10 and 41. 
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The upper limit for the potential rain pH ordinate is based upon an assumed 
rainfall rate of 25 mm/hr. The areas beyond X = TOO km downwind and 10 hr 
postlaunch should be treated with an additional level of uncertainity, as dis- 
cussed previously. Thus the respective inclusive envelopes of these results 
define a regime of most probable initial rain (potential) pH when it is recog- 
nized that the model application is restricted to conditions of low-to- 
modecate SRM cloud humidity (HCA mostly in gas phase) and to ground locations 
beneath the SRM CC on which an independently generated overriding rainfall 
commences under stable stratification conditions aloft. 


Examples of Acidic Rain Deposition Characteristics 

Before detailed examples of calculated rain pH and HCil deposition due to 
SRM cloud washout are presented, consider briefly the characteristics of natural 
rain and contrast these with the regime of concern for individual launch 
events, where potential rain pH may be <3 and sometimes <1. Figure 18 shows 
some characteristics of rain acidity for the Florida area (refs. 44 and 45) and 
also a set of parametric relationships among potential pH, effective SRM cloud 
thickness, and vertically averaged HCJl(g) concentration. The pH was calculated 
from equation (46) at a conservative rainfall intensity of 25 mm/hr. 

The volume-weighted average rain pH in the Cape Canaveral area was deter- 
mined to be 4.61 ± 0.22 (monthly standard deviation), based on an extensive 
study of rainwater characteristics (ref. 45) over the period July 1977 to June 
1979. Although the range of rain pH for north-central Florida during the 1967- 
1968 period (5.3 to 6.8) is based on only four measurements (ref. 44) and is 
therefore not necessarily trepresentative, these unusually high values are con- 
sistent with a downward temporal trend in deduced for the northeastern 
United States over the last 25 years (refs. 46 to 50) . More specifically, a 
recent comprehensive study of rainwater characteristics in Florida (ref. 51) 
leads to the conclusion that a significant decline in pH and bicarbonate con- 
centration has occurred for northern and central Florida over the past two 
decades. 

The open circles in figure 18 represent potential pH corresponding to the 
source- strength- modified SFW Shuttle case illustrated in figure 15. Obvi- 
ously, these predictions of single-event rain acidity are significant even for 
relatively large distances when compared with the pH of natural precipitation. 

Detailed predictions of downwind pH (eg. (50) or (51)) and cumulative HC?, 
deposition (eg. (58) or (59)) are illustrated in figures 1 9 to 24 for 

the FFW Titan III, SFW Titan III, and fully modified SFW Shuttle cases. Poten- 
tial quantities, pHpot fi’oni equation (52) and Gx,pot equation (61), 

which apply at the first onset of rain are shown as dashed lines. Examples of 
progressive washout for various assumed onsets of rainfall at each of three 
different intensities are shown as solid lines. Cumulative areal deposition 
applies along the CC path; off-axis deposition can be calculated from equa- 
tion (60) , with the effective SRM cloud diameter D^ evaluated from equa- 
tion (55). Note that attainment of CO 2 equilibration at pH = 5.7 implies 
that 100 percent washout of the HCi has occurred upwind. Nearly complete wash- 
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out can readily occur within 30 km at high rainfall rates. Finally, the corre- 
sponding downwind growth of SRM cloud diameter is shown in figure 25 for each 
of these three cases. 

Several observations of the pH and results shown in figures 19 to 24 

are now cited briefly. While rain acidity decreases with increasing rainfall 
rate, total deposition increases. Also, high rainfall rates occurring close to 
the launch site result in significantly greater depositions than lower rainfall 
rates occurring farther out. The progressive increase in downwind potential 
rain pH, due exclusively to dispersive decay of HCA column density for any given 
rainfall rate, is much slower for the FFW case than the SFW Titan III case. 

These cases constitute the two limiting dispersion extremes with respect to 
downwind distance. Note that for both the FFW and SFW cases the pH of a light 
drizzle (1 mm/hr), beginning at 15 km downwind and continuing out to 100 km, 
increases very slowly with respect to potential pH (e.g., 0.6 to 0.8 pH units), 
and is about equivalent to the potential pH for heavy rainfall (25 mm/hr) onset 
at «80 km. 

While the main thrust of this section has centered on the prediction of 
acid rain characteristics, some closing remarks on the possible consequences of 
SRM cloud washout are now given. In general, rain acidity derived from washout 
of an overhead SRM cloud will tend to determine the respective transient rates 
of acid reaction with various exposed surfaces, either manmade or natural. The 
duration (and acidity) of any primary exposure, and hence degree of reaction 
during exposure, will depend on SRM cloud characteristics, wind speed, and loca- 
tion. Subsequent cleansing of nontrapping surfaces and dilution of accumulated 
acid rain will occur if rainfall continues after SRM cloud passage. In this 
case quenching or reduction of the rate of acid attack will tend to occur, the 
extent depending partly on the nature of each surface. However, regardless of 
any subsequent cleansing phenomena, and also recognizing the low probability of 
event reoccurrence, one must still consider the possibility that localized 
exposures to rain containing HCil and aluminum oxychlorides may have both short 
and long term effects on the local ecology (refs. 39 and 40) . Long term damage 
effects may depend not only on the concentration history of acidic rain, but 
also on the deposition histories of both acidic rain and subsequent rainwater 
dilutions, as well as total accumulated deposition and a host of other environ- 
mental parameters. Further guidance in defining possible long term effects may 
be obtained from recent studies of regional acid rain problems in parts of 
Europe and the northeastern United States (refs. 46 to 51). Although some of 
the more obvious and serious environmental effects have already been defined 
in these studies, it is clear that understanding of the problen is far from 
complete. 


Dispersive Decays of Average HCA Concentration 

Dispersive decays of corresponding vertically averaged HCil concentration 
were deduced for the eight Titan III cases, and the results are shown as func- 
tions of X and elapsed time in figures 26 and 27, respectively. These 
averaged peak-center-line UCH concentrations were computed from expressions 
for o by (a) estimating effective cloud thicknesses which applied for 
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Xcs = 5 km as indicated in table IV, and (b) dividing the appropriate expres- 
sions for o in table IV by the corresponding cloud thicknesses. 


Since the dispersive decays of ECU concentration in figures 26 and 27 cor- 
respond to decays of potential rain pH in figures 16 and 17, comparisons of 
these figure sets facilitate preliminary assessments of acid rain potential 
based on measurements of in-cloud HCil. Such comparisons have proven useful in 
correlating in-cloud HCS. concentrations with simultaneous rain composition data. 
(See ref. 43.) 

A comparison of the calculated results in figure 27 with experimental mea- 
surements of peak in-cloud HCA concentrations for eight Titan III launches is 
illustrated in figure 28 (reproduced from ref. 43). Since the calculated dis- 
persion results were derived from a set of standard meteorologies for the Cape 
Canaveral area, they have no direct relationships to the actual launch cases. 
However, certain first-order features of the respective sets are in agreement. 

First, to a good approximation, all the experimental in-cloud HCA concen- 
tration data appeared to be adequately characterized by single-term power-law- 
decay expressions, such as shown to apply for the modeled results. Next, the 
total range of measured HC£ concentrations is approximately bounded by the 
total range of predicted concentrations after «0.2 hr (lower limit of model 
validity) and up to the indicated termination of each launch data fit (at 0.5 
to 5 hr). Moreover, the respective ranges of highest to lowest HCil concentra- 
tion were large; e.g., the ratios of highest to lowest HC?, were «100 for both 
sets and increasing after 1 hr postlaunch. Despite this intersection of the 
experimental and calculated sets of HCil decay histories, four of the measured 
in-cloud peak HCil concentrations decayed significantly slower with time than 
the slowest calculated (g = 0.83 for FFW) . Notably, the worst-case model cal- 
culations of HCil concentration (highest values) were partly confirmed experi- 
mentally. The apparent tendency of the dispersion model to overpredict the 
rate of decay of peak HCil concentrations is consistent with the derivation of 
MDM empirical dispersion coefficients, which are essentially based on correla- 
tions of relatively small-scale turbulence measurements. Thus, while there is 
surprisingly good overall agreement with respect to HCil concentrations, the 
present MDM application is believed to be deficient in that it fails to account 
for the reduced effect of large-scale turbulent diffusion, bounded by an inver- 
sion layer, that inevitably applies to large stabilized SRM exhaust clouds sev- 
eral minutes after launch. 


Test of Precipitation Scavenging Subroutine in MD^^5(II) 

The present analytic HC)l(g) washout model allows an independent check on 
the output of the multilayer (and essentially numerical) calculations of 
MDM-5(II). The comparison made here was to examine the progressive washout of 
HCil from an SRM cloud, as now described. (Deposition calculations were not 
available.) The FFW meteorology was selected; a washout coefficient correspond- 
ing to 7.7 mm rain/hr was specified (A = 4.68 x 10“^ sec”^); and the elapsed 
time from cloud stabilization for onset of rain was defined (t*| = 394 sec) . Two 
calculations of vertical HCil(g) profiles were needed, one with MDM-5(II) that 
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allowed evaluation of a (with rain) and a reference calculation with MDM-4(II) 
that defined a (no rain). The latter was equivalent to the results presented 
in figure 16. Equating a (with rain)/a(no rain) with the ratio p(HC2/)/Pq(HC 2/) , 
and plotting selected results as a function of distance fron cloud stabilization 
X^g led to the results shown in figure 29. The straight line, which agrees 
very well with the calculated points, was obtained from an alternate form of 
equations (49), i. e. , 


Note that p(HC8/)/Pq (HCH) =1 at = 2.41 km, which corresponds to the onset 

of rain at t*] = 394 sec. Results of a similar test of MDM-5(II) using the SFW 
meteorology are shown in figure 30 and indicate equally good agreement. The 
difference in downwind HCil(g) washout histories between figures 29 and 30 is due 
entirely to the difference in average SRM cloud speed U^. 

These comparisons demonstrate that calculated HCA(g) washout, evaluated by 
equations (9) and (49) for the present analytic model, was quantitatively equiv- 
alent to that evaluated numerically by using MDM-5(II). As a final note it is 
suggested that the simplicity of the present model, which requires only a single 
determination of a decay in the absence of rain, allows subsequent evalua- 
tions of a and various HCA washout and deposition characteristics with greatly 
reduced computational effort. Thus, a comprehensive matrix of lengthy MDM-5(II) 
calculations is no longer required to evaluate the effects of various combina- 
tions of assumed rain onset time and washout coefficient corresponding to 
assumed rainfall intensity. It is conceivable that this might prove valuable 
for real-time predictions and probability assessments of acid rain hazards dur- 
ing prelaunch countdowns. 


p(HCil)/Po(HCil) = exp 


-A (1000 t-| 

V Uc 


SUMMARY OF RESULTS 

A previously developed HCS/(g) washout model for rain scavenging of solid 
rocket (SRM) exhaust clouds was refined and applied to nine independently 
determined SRM exhaust cloud dispersion cases. The model refinements consist of 
(a) including more representative values for air properties and diffusion 
coefficient; (b) developing a more comprehensive analytic approach for general- 
izing HCi^(g) absorption, washout, and resultant rain pH and HCi^ deposition char- 
acteristics as functions of droplet diameter, terminal velocity, raindrop size 
distribution and intensity, and exhaust cloud dispersion characteristics; 

(c) deriving two new HCil(g) washout coefficient es^ressions, based on an experi- 
mentally verified convective diffusion equation (modified Frossling) and inte- 
grated over both the Marshall-Palmer (M-P) raindrop size distribution and the 
previously used Kelkar raindrop size- intensity data; and (d) deriving a new 
HCil(g + aq) washout coefficient es^ression for SRM exhaust aerosol, based on a 
detailed review and analysis of published laboratory data on rain scavenging of 
SRM exhaust and integrated over the M-P distribution. 
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The model treats the idealized case of an independently generated vertical 
rainfall that overrides and scavenges an independently dispersing SRM exhaust 
cloud under stable stratification conditions in the lower troposphere. The 
first two relationships for washout coefficient A versus rainfall intensity 
H derived in this paper characterize washout of HCS,(g). These apply at low- to- 
moderate ambient relative humidities, where it can be assumed that HCJl(g) will 
predominate over the acjueous acid aerosol component after a few minutes of SRM 
cloud dilution. The geometric mean of these two power-law expressions, 

A = 1 . 39 X sec“^, is recommended for subsequent predictions of HCil{g) 

washout under idealized average conditions. 

The third A expression characterizes HCil (g + aq) washout for diluted 
and humidified SRM exhaust aerosol. This expression has the same analytic form 
and is both statistically and numerically equivalent, to a good approximation, 
to that obtained for HCil(g) washout. Thus it is provisionally recommended that 
the geometric mean HCil(g) washout coefficient should also be applied to charac- 
terize HCA(g + aq) washout from SRM exhaust clouds under moderate ambient 
humidity conditions {^95 percent), where significant quantities of HCJl (aq) may 
coexist with HCil(g). 

The washout model has been generalized for application to SRM exhaust 
clouds as follows. The resultant p*[ of rain and HCi?/ deposition rate are char- 
acterized in terms of (a) the initial source strength of HCil; (b) a power-law 
decay of vertical HCil(g) column density a = aX"^, where X is downwind dis- 
tance from the launch site and a and 3 are empirically determined constants; 
and (c) the downwind distance (time) at which a steady overriding rainfall com- 
mences at a specified intensity. The cumulative deposition of HCil is then char- 
acterized in terms of ground coordinates when (a) to (c) are specified and the 
SRM cloud geometry is assumed to be an expanding right circular cylinder. 

The washout model was applied to a refined "spring fair weather" (SFW) 

Space Shuttle exhaust cloud dispersion case and eight Titan III dispersion cases. 
All nine cases were based on independently calculated downwind dispersive decays 
of a, deduced previously fron application of the NASA Marshall Space Plight 
Center Gaussian multilayer diffusion model, MDM-4(II), to seven standard meteo- 
rological regimes. These regimes were originally selected to represent a range 
of characteristic cases of overland advection with turbulent diffusion in the 
planetary boundary layer for the Cape Canaveral, Florida, area. 

Results obtained from application of the washout model to the eight 
Titan III cases indicated that the dispersive decays of potential rain pH, i.e. , 
volumetric average pfi at the first onset of rainfall, differed greatly among 
the seven standard meteorological regimes. A range of more than 2 pH units was 
spanned at X ^ 1 00 km downwind and/or t ^ 2 hr postlaunch elapsed time. At 
shorter distances and times, the total span was less but still exceeded an order 
of magnitude for X > 1 0 km and t > 0. 2 hr. These large spreads in potential 
pR are consistent with earlier expectations of atmospheric dispersion under 
widely different meteorological conditions. They also highlight the need for 
developing a more realistic and comprehensive atmospheric dynamics (with cloud 
microphysics-scavenging processes) model to deal effectively with this uncer- 
tainty. Environmentally significant potential rain acidities (pH ^ 1.5) for 
one-time exposures were shown to be possible (for H ^ 25 mirv'hr) at downwind 
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distances up to 50 km and elapsed times of 5 hr for the two most severe, least 
dispersive Titan III cases. This result contrasts with a measured volumetric 
average pH of 4.61 ± 0.22 (monthly standard deviation) for naturally occurring 
acidic rains in the Cape Canaveral area over the last 2 years. 

Detailed examples showing calculated downwind histories of rainwater pH. and 
HCil deposition (for both potential and progressive washout) were shown for the 
least and most dispersive Titan III cases (fall fair weather (FFW) and SFW) and 
a modified SFW Shuttle case. In general, rain acidity decreased with increas- 
ing rainfall rate while deposition increased; this is a consequence of the <1.0 
power-law dependence of washout coefficient on rainfall rate. High rainfall 
rates occurring close to the launch site resulted in much greater depositions 
than lower rainfall rates occurring farther out. Nearly complete washout of HC^^ 
from an SRM cloud could occur within 30 km at high rain intensities. Progressive 
washout at a low rainfall intensity (e.g., 1 rm/hr) over large distances (e.g., 

15 to 100 km) results in relatively slow decay of rainwater acidity (0.6 to 
0.8 pH unit). Thus, damage to ground-receiver surfaces may tend to be greater 
(Icwer pH) and more extensive (longer footprint) at lower rainfall intensities. 


CONCLUDING REMARKS 

Several factors affect the validity of the predicted £51 values and ground 
depositions. The greatest uncertainties stem from the basic idealized assump- 
tions used to formulate the probleu. The assumptions of an independently 
dispersing SRM cloud that experiences washout from an independently generated 
overriding rainfall under stable stratification conditions specifically exclude 
(a) significant convective activity associated with storm dynamics and (b) rain- 
out processes stemming from the interaction of SRM cloud components, HCS/ (g) and 
HCA (aq) on chi orided- alumina nuclei, with natural clouds. Since rainfall in the 
Cape Canaveral area is frequently associated with significant convective activ- 
ity, it seems imperative that a more realistic atmospheric dynamics model be 
developed with appropriate par ameterizat ions of the essential cloud microphysics- 
scavenging processes. 

Additionally, uncertainties related to the input o decays also affect the 
validity of the predictions, especially at large distances (e.g., X ^ 100 km). 

A refined inclusion of (a) convective loss of HCit from the SRM cloud's upper 
boundary, (b) HCA sorption at ground level, and, probably, (c) a more realistic 
treatment of horizontal wind shear effects for each layer would have increased 
potential rain Since the decays of cr stemmed from MDM-4(II) calculations 

based on relatively small-scale turbulence measurements, use of variances scaled 
more appropriately to cloud size would have resulted in smaller but more realis- 
tic values of 3 at large X and t, similar to those observed for Titan III 
clouds. The effect in this case would be to decrease potential rain pH at 
large X. Clearly these competing effects require complex analyses beyond the 
scope of this paper. 

Finally, the rain scavenging relationships derived in this paper apply for 
washout of HCil(g) and HCA(g + aq) by rain of specified intensity and average 
polydispersity. Naturally occurring rains vary widely in both respects, and 
their characteristics are very difficult to predict. Rainout of HCJt associated 
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with convective cloud activity may be important or dominant under many condi- 
tions, but these processes have not been examined herein. Although HCA (g) 
will tend to predominate over HCH (aq) on chlorided-alumina nuclei after a few 
minutes of cloud dilution at low to moderate ambient humidities, the applica- 
bility of the deduced HCA (g + aq) washout coefficient is most questionable at 
very high relative humdities, say ^95 percent, which frequently exist at Cape 
Canaveral. Thus since washout coefficients for aerosols depend uniquely and 
significantly on aerosol size distribution and since the distribution of aero- 
sol sizes and the HCS- (aq) concentration exhibit complex microphysical time 
and concentration dependencies, refinements in the description of HCJl (g + aq) 
washout may require detailed cloud microphysics calculations. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 13, 1981 
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REVIEW OF SCAVENGING STUDIES 
Knutson and Fenton 

In the study of Knutson and Fenton (ref. 24) 0,09-cm water drops, at 
terminal velocity, scavenged the exhaust of a small, solid-propellant rocket 
motor in a 5 . 5-m-diameter spherical chamber. Bubblers were used to determine 
chamber HCil concentration, which included HCA(g) and HCiL{aq) aerosol attached to 
^^2^3 particulates. A log- log plot of CA* scavenged against chamber C4“ concen- 
tration indicated approximately a unity slope based on eight data points. How- 
ever, the resultant proportionality constant was approximately 0.5 times the 
value predicted by the modified Frossling correlation (eq. (3)). These authors 
also observed that relative humidity in the range 69 to 98 percent had no effect 
on scavenging rate, based on this plot. Relative humidity was calculated from 
total initial water content, at measured dry bulb temperature, without account- 
ing for subsequent HCi (aq) aerosol formation and solute vapor pressure lowering 
effects. The authors calculated a washout coefficient expression based on drop- 
let HCil(g) absorption at 0.5 times the rate expressed by the modified Frossling 
correlation, the M-P raindrop size distribution, and = 0.173 cm^/sec 

at 15^ and 1 atm. Numerical integration and empirical fitting of the results 
led to their recommended expression for HCt washout (ref. 24); 


A = 8.3 X 10"5h0*567 


(Al) 


Unfortunately, the significance and utility of the Knutson and Fenton wash- 
out coefficient (eq. (Al ) ) is unclear. First, the effective air-exhaust mass 
ratio was about 225 for six of the eight useful data runs. This contrasts with 
exhaust cloud ratios of ^10^ after a few minutes coupled with a much greater 
availability of water vapor in the real case. Thus, co-condensation of HC8/(g) 
and H20(g) in the chamber experiments, forming HCMaq) on chlor ided-alumina 
nuclei (see refs. 3, 52, 53, and 54) and on the walls (discussed in ref. 24) , 
favored large HC4 partitioning and water vapor pressure reduction effects. 

That is, a relatively high concentration of hygroscopic particles competed for 
a relatively limited supply of water vapor. 

As a check, the present author calculated the properties of a flat-surface 
(no Kelvin effect) locally equilibrated HCil(aq) aerosol based on the cited 
experimental data and an iterative analytical solution of the respective parti- 
tioning equations for HCil and H 2 O by using the methodology outlined in refer- 
ence 43. For run 7 in reference 24, equilibrium HCA(aq) aerosol formation would 
reduce the relative humidity from 98 to 86.4 percent and would leave only 
7.76 ppmv HCil(g) out of the measured 200 ppnv HCil(g + aq) during the rain scav- 
enging period. Similarly, for run 8 in reference 24, relative humidity would be 
reduced frcxn 87 to 83 percent, thus leaving only 16.5 pfanv HCft(g) out of the 
measured 100 ppmv HCA(g + aq) . The respective molalities of HCi(aq) are 3.08 
and 3.64 for runs 7 and 8. These results indicate that only small percentages 
of the chamber HCjL existed as HCa(g) (3.9 and 16.5 percent, respectively) and 
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the remainder existed as concentrated Hdl (aq) aerosol with a relatively fine 
but essentially unknown size distribution. Therefore, since washout of Hdl (aq) 
aerosol by raindrops is typically an order of magnitude less efficient for 
micrometer- sized droplets than for Hdl (g) but is strongly dependent on aerosol 
size distribution (refs. 14 and 17), the data of reference 24 are of limited 
usefulness, even though the observed overall trend (0.5 times A for Hdl (g) ) 
is in the e3q>ected direction for moderately small (e.g., 1- to 5-)im diameter) 
aerosol sizes. Parenthetically, note that a predominance of submicrometer aero- 
sol sizes could decrease this factor to values approaching 0.1, while a predomi- 
nance of lO-jim aerosol sizes could increase the factor to »1. (See ref. 14.) 


Fenton and Purcell 

The follow-on IITRI study by Fenton and Purcell (ref. 25) was an attempt 
to overcome some of the experimental difficulties of the first study; a number 
of ancillary measurements were also conducted. In particular, they sought to 
conduct rain scavenging measurements on laboratory- produced rocket exhaust 
clouds at lower Hdl concentrations than before, typical of poststabilized SRM 
clouds. The solid propellant used was essentially the same; the rocket motor 
and propellant mass were similar but about 50 percent smaller; and the nozzle 
was identical. 

The rocket motors were fired into the same spherical (5.5 m) test chamber 
used previously, but the e3q>er imental approach differed significantly from this 
point on. At preselected times partly diluted SRM exhaust in the primary 
sjAierical chamber was conducted into a Teflon^ 1.5-m^ experimental chamber via 
a specially designed aerosol transport tube, which incorporated a peripheral 
transpiration air sheath and allowed further dilution and efficient transport 
of the aerosol. 

The raindrop generator used to produce uniform-size droplets was essen- 
tially the same as before (ref. 24). However, three different droplet sizes 
(0.055, O.n, and 0.30 cm) were utilized instead of the single size (0.09 cm) 
used previously. Terminal velocity was claimed to have been achieved at the 
raindrop generator exit for the two smallest sizes but not for the 0.30-cm 
droplets. A glass funnel attached to an N 2 purged holding flask was used to 
collect the rain. This replaced the mineral oil collection pan used previ- 
ously. Exit ports located along the periphery of the funnel were used to 
remove a matched flow of N 2 sheath gas. Parenthetically, a denser gas such 
as argon would have been superior. 

The concentration of HCil in the test chamber was measured with a Geomet, 
Inc., chemiluminescent detector rather than an aqueous bubbler (impinger) system 
used previously. A single comparison of the two methods indicated that the bub- 
bler system measured about 10 percent more chloride. This small difference may 
not have been experimentally significant, but it was consistent with expecta- 
tions of additional aluminum oxychloride salts trapped as chlorided-alumina par- 
ticulates by the bubbler. This result coupled with Geomet HCS. (aq) calibration 


^Teflon: Registered trademark of E. I. du Pont de Nemours & Co., Inc. 
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checks appeared to confirm that the Geomets used to monitor HC£ in the Teflon 
test chamber were measuring HCil (g + aq) but not chloride associated with the 
alumina. 

The results of 33 rain scavenging tests are summarized in table VI of ref- 
erence 25, As a check the present author independently calculated values 

by (a) assuming various effective scavenging heights consistent with the afore- 
mentioned setup, (b) using the tabulated rain collection data to evaluate cor- 
rected raindrop chloride concentrations to at least two significant figures, and 
(c) using the cited average raindrop diameters. An assumed scavenging height 
of 145 cm was found to give reasonable agreement with all but four of the tabu- 
lated values (2a, 2b, 4a, and 4f in table VI of ref. 25) • These appeared 

to represent calculation- typographical errors, and with one exception (4f), were 
not obvious misfits in their final parameterized-data correlation (fig. 18 of 
ref. 25). 

Fenton and Purcell sought to generalize these scavenging data by defining 
a quasi-dimensionless group of variables, ^A/Pair^c ^d2, in g HCH absorbed/g 
air. They plotted (on a log-log scale) this parameter as a function of drift- 
corrected HCil (g + aq) concentration, C(HCil). Note that in the absence of a 
more specific definition for C(HCil) it must be assumed that pprav units apply 
to C(HCil) . 

Although the parameterized-data scatter was relatively large in figure 18 
of reference 25, best-fit empirical expressions were stated to be 


^A 

= 2.82 X 10“7c(HC!1)0*®24 

PairV<„d2 

= 6.72 X 


For the latter, was defined as the equivalent mass concentration of HC^l, in 

g HC£/m^. Clearly the first expression, in terms of C(CHJl), does not agree 
(factor of «1.2 to «2.5 high) with the best-fit line drawn in figure 18 of ref- 
erence 25. Additionally, it can be shown that the latter expression in terms of 
is inconsistent with the first expression when an appropriate value for lab- 
oratory air density is used and is in greater disagreement with the best-fit 
line shown in figure 18 of reference 25 than the C(CHil) expression. Reasons for 
these disagreements are not apparent. 

Fenton and Purcell then evaluated A as follows. They used the expression 
for parameterized raindrop scavenging rate to evaluate a washout coefficient foe 
polydisperse rain. They defined an integral equation equivalent to 


* If <3(n) 

A = — \ Ra d(d) (A2) 

Ma J d(d) 
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^ ^ ■'0-5MA°*®2^)(PairV„d2)ii^ d (d) (A3) 


Fenton and Purcell used the M-P raindrop distribution (eq* (25)) and an approx- 
imate power-law expression for droplet terminal velocity to perform numerical 
integrations for calculation of A as a function of H. These values of A 
were stated to be best fit by the expression (incorrect, exponent for 
should be negative) 


A = 4.21 X 


Using Pair “ 1216 g air/m^ for their standard conditions, they obtained a 
final expression (also incorrect) for the washout coefficient; 
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STANDARD METEOROLOGICAL REGIMES FOR CAPE CANAVERAL 

The meteorological regimes for this work consist of the entire set of stan- 
dard meteorologies (ref. 5) that were originally used to assess atmospheric dis- 
persion of rocket exhaust effluents in the Cape Canaveral, Florida, area. These 
regimes represent the major meteorological conditions which lead to overland 
transport and are likely to be encountered in the Cape area. 

The meteorological profiles of temperature, wind speed, and wind direction 
for all seven meteorological regimes are shown in figures 11(a) to 11(g), as 
reproduced frcm references 33 and 28 (SB) . The profiles are based on averages of 
selected data obtained from Cape Canaveral rawinsonde releases and from the NASA 
150-Meter Meteorological Tower at Kenndey Space Center, as tabulated in refer- 
ence 5. Inspection of these profiles indicates large differences in the thermal 
structure and wind characteristics aloft. Note that significant deviations 
occur, at various levels, fran the preprocessor-calculated average values of 
wind speed and wind direction 0^ (summarized in table III of main body of 

this paper), that were used in MDM-4(II) to describe average SRM cloud advective 
motion after altitude stabilization. Thus, horizontal wind shear effects are 
not directly accounted for in the MDM-4(II) calculations since individual layers 
in the multilayer structure are not allowed to advect independently, even though 
average values for each layer are specified in the program output. Instead, the 
effects of horizontal and vertical wind shear, relative to average SRM cloud 
advective speed and direction, are simulated in each layer by respective sets of 
horizontal (cross direction) and vertical turbulence parameters applied to each 
layer. These parameters are standard deviations of the wind azimuth 
and elevation angle Og fluctuations. While these standard deviations can be 
obtained from a comprehensive set of local meteorological data for the Cape 
Canaveral area (refs. 34 to 36), they can also be deduced from the profile mea- 
surements, such as those shown in figures 11(a) to 11(g). The latter procedure 
was used for the seven standard meteorologies employed in this study; the meth- 
odology is described in references 26 to 32 and the specific parameters used are 
tabulated in an unpublished paper by G. L. Pellett. 
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CHARACTERISTICS OF ALTITUDE-STABILIZED SRM SOURCE CLOUDS 

An example of the geometric definition of a Titan III stabilized SRM 
exhaust cloud is shown in figure 12(a) for the sea breeze (SB) meteorological 
regime at Cape Canaveral, as obtained frcxn reference 28. The calculated height 
of the cloud centroid Zj^ = 832 m is based on an early value for the effec- 
tive heat release Qi = 691 cal/g propellent; 2790 cal/g is used presently. 

The fine-layered surface mixing layer has an upper bound (Hj^ = 800 m) based 
on the onset of temperature inversion at 800 m as shown in figure 11(d). The 
stem of the so-called column cloud begins at 1300 m with a radius of 200 m. 

The reason for assigning 1300 m as the transition from the so-called stabi- 
lized ground cloud to the stabilized column cloud is apparent in figure 11(d), 
where the temperature inversion terminates and the gradient of local (and 
potential) temperature turns negative. 

A second stabilized SRM cloud is shown in figure 12(b) for the post-cold 
front passage meteorological regime applied to the case of an on- pad abort of 
the Titan III vehicle (Post-CF, pad abort). In this case, obtained frcxn refer- 
ence 28 and based on Qj = 691 cal/g, the calculated height of the CC is 
Zjn = 1132 m, compared with 751 m for the normal launch case (Post-CF) and 832 m 
for the previously discussed SB case. Thus, both the vertical temperature pro- 
file and the input source strength influence the determination of Zj^. Note 
also that the height of the surface mixing layer (Hj^^ = 1400 m, fig. 11(g)) is 
greater than in the SB case but is not defined strictly at the minimum tempera- 
ture point (1700 m) in figure 11(g). Instead, Hj^ is defined at the turning 
point in potential temperature. 

A third Titan III SRM cloud, also obtained from reference 28 and based on 
Qj = 691 cal/g, is shown in figure 12(c) for the cold front passage (CFP) case. 
The corresponding vertical temperature gradient in figure 11(f) was not so 
steep as in some cases, resulting in a relatively low CC height (Z^ = 675 m) . 
However, the temperature profile did not exhibit an inversion up to 2000 m, and 
thus the surface layer was assigned a depth of Hjj^ = 2000 m; the multilayer 
structure was arbitrarily terminated at this point. Geometric definitions for 
the remaining meteorological cases were based on application of the same prin- 
ciples first described. Plots of the layer structures and initial vertical 
distributions of HCS^ for the present cases (Qj = 2790 cal/g propellant) are 
similar, based on this author's unpublished tabulated data. 
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Figure 1 Empirical fits of Sh/d^V^ parameter for convective diffusion 
of HCil(g) to raindrops (modified Frossling) falling at terminal velocity 
in air at 1 5^C and 0.85 atm. 
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Figure 2.- Empirical fit of d(Sh) parameter for convective diffusion of 
HCii(g) to raindrops (modified Frossling) falling at terminal velocity 
in air at 1 5^C and 0.85 atm. 
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Figure 3.- Empirical fits of HCil(g) washout coefficients based on modified 
Frossling equation and applied to Marshall-Palmer (M-P) raindrop size 
distribution function and Kelkar's raindrop intensity data. Expressions 
apply at 1 5^C and 0.85 atm. 
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Figure 4.- Empirical fit of weighted mean Sh/d^V^ parameter for convective 
diffusion of HCil(g) to polydisperse raindrops (modified Frossling) falling 
with intensities measured by Kelkar for average rains. Expression applies 
at 1 5°C and 0.85 atm. 
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Figure 5*- Empirically deduced relationships for HCA (g + aq) washout based 
on Fenton and Purcell's (F-P) data (ref. 25) on scavenging of diluted SRM 
exhaust in laboratory chamber by uniform-size raindrops. They include 
power-corrected F-P expression (eq. (35)) and present independently 
deduced expression (eq. (40)). 
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Figure 6,- Analysis of Fenton and Purcell's data (ref. 25) on laboratory 
scavenging of diluted SRM exhaust by 0 . 055-cm-diameter water droplets 
falling at terminal velocity. Solid line represents mean value of 
Ra/C(HCA). Numbers represent relative humidity. 



Droplet scavenging rate for HC£, R. , g HC£/sec 


□ Data for 0.11-cm droplets 


58 / 
sy D58 
^□66 
□ 66 
□ 60 


Concentration of HCf (g + aq) 
in chamber, C(HCj^), ppmv 


Figure 7.- Analysis of Fenton and Purcell's data (ref. 25) on laboratory 
scavenging of diluted SRM exhaust by 0.11 -cm-diameter water droplets 
falling at terminal velocity. Solid line represents mean value of 
R;^/C(HCil). Numbers represent relative humidity. Data points with 
asterisks had relatively large chloride corrections and are considered 
less reliable. 
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Figure 8.- Analysis of Fenton and Purcell's data (ref. 25) on laboratory 
scavenging of diluted SRM exhaust by 0. 30-cm-diameter water droplets 
falling at accelerating subterminal velocity. Solid line represents 
mean value of Ra/C(HCA). Numbers represent relative humidity. Data 
points with asterisks had relatively large chloride corrections and 
are considered less reliable. 
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Figure 9.- Effect of droplet diameter on concentration-normalized HCit (g + aq) 
scavenging rate on a per-droplet basis. Least-squares fit of log [r^/C(H 
versus log d is shown for mean values of Rj^/C(HCit) illustrated in fig- 
ures 6 to 8. Vertical bars correspond to standard errors of mean values. 
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Figure 10.- Empirical fit of HCA(g + aq) washout coefficients (circles), based 
on present analysis of Fenton and Purcell's SRM exhaust rain scavenging 
data (ref. 25) and Marshall-Palmer (M-P) raindrop size distribution. Equa- 
tion (26) for HCS-(g) washout, based on modified Frossling equation and 
M-P distribution, is shown for comparison. 
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Figure 
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For fall fair weather (FFW) meteorological regime (from ref. 33) . 

n Vertical profiles of temperature, wind speed, and wind direction 
at Cape Canaveral, Florida. 
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(c) For low-level sea breeze (LLSB) meteorological regime (from ref. 33). 


Figure 11.- Continued. 
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(d) For sea breeze (SB) meteorological regime (from ref. 28) . 


Figure 11 Continued. 




Height, 



(e) For fair weather, pre-cold front (FW, Pre-CF) meteorological regime 

(from ref. 33) . 


Figure 11.- Continued. 
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Figure 11 Concluded. 
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(a) For sea breeze (SB) meteorological regime. Height of cloud centroid is 
832 m, and surface mixing- layer depth is 800 m (from ref. 28) . 

Figure 12.- Dimensions of stabilized ground cloud of Titan III exhaust products 

(Qj = 691 cal/g) at Cape Canaveral, Florida. 
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(b) For post-cold front meteorological regime during normal launch (Post-CF) 
and during on-pad abort (Post-CF^ pad abort) (from ref. 28) , 


Figure 1 2.- Continued. 
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Figure 13.- Sane HCi^- source-strength characteristics for Space Shuttle and 
Titan III SRM launches. Note that cumulative source strength in MDM-4(II) 
stabilized cloud is less than cumulative line source input. 
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Figure 14.- Functional relationship used in reference 3 to characterize 
downwind dispersive decay of a for Space Shuttle SRM exhaust cloud 
under SFW meteorological conditions. HCil source strength was 
iDo = 89 X 10^ g HCil. 
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Figure 15.- Virtual-source-distance corrected empirical functions which 

characterize modified downwind dispersive decay of a for Space Shuttle 
SRM exhaust cloud under SFW meteorological conditions. Functions for 
both original (89 x l 0^ g HCA) and presently modified (61 x 1 0^ g HCS.) 
stabilized-cloud source strengths are shown. 
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Lgure 16.- Dispersive decay of a a 
distance, for eight Titan III case 
{dashed line) of exhaust cloud dis 
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Figure 17.- Dispersive decay of CJ and potential rain pH as function of downwind 
drift time, for seven Titan III cases (solid lines) and one Space Shuttle case 
(dashed line) of exhaust cloud dispersion. Results are based on application 
of MDM-4(II) to seven standard meteorological regimes, equation (52) for 
and equation (42) for A (gectnetric mean) . 
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Figure 18.- Characteristic rain acidity for Florida area and relationships 
among potential pH, effective SRM cloud thickness, and vertically averaged 
HCH(g) concentration. The latter case was calculated for rainfall rate 
of 25 mm/hr using equation (46) with t = 0. 
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Figure 19.- Calculated volumetric-average rain acidity due to HCJ!-(g) washout for 
FFW Titan III case. Potential rain pH (dashed line) applies at first onset 
of rain; pH of steady rain (solid line) reflects progressive washout, which 
was begun at 15, 30, and 50 km downwind. 
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Figure 20.- Calculated cumulative HCil deposition per unit area along cloud- 
centroid path due to HCil(g) washout for FEW Titan III case. Potential 
deposition (dashed line) applies at first onset of rain for entire SRM 
cloud. Deposition after start of steady rain (solid line) reflects 
progressive washout, which was begun at 1 5, 30, and 50 km downwind. 





<u 

> 

(D 

d 

O 

u 

b/3 


0) 

-i-j 

C 

C4H 

o 

a 

a 



Figure 21.- Calculated volumetric-average rain acidity due to HCil(g) washout 
for SFW Titan III case. Potential rain pH (dashed line) applies at first 
onset of rain; pH of steady rain (solid line) reflects progressive washout, 
which was begun at 15, 30, and 50 km downwind. 
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Figure 22.- Calculated cumulative HCil deposition per unit area along cloud- 
centroid path due to HCil(g) washout for SFW Titan III case. Potential 
deposition (dashed line) applies at first onset of rain for entire SRM 
cloud. Deposition after start of steady rain (solid line) reflects 
progressive washout, which was begun at 15, 30, and 50 km downwind. 
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Figure 23.- Calculated volumetric-average rain acidity due to HCil(g) washout 
for source-strength-modified SFW Shuttle case. Potential rain pH (dashed 
line) applies at first onset of rain; pH of steady rain (solid line) reflects 
progressive washout, which was begun at 15, 30, and 50 km downwind. 
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Figure 24.- Calculated cumulative Hdt deposition per unit area along cloud- 
centroid path due to HCit(g) washout for source-strength-modified SFW 
Shuttle case. Potential deposition (dashed line) applies at first onset 
of rain for entire SRM cloud; deposition after start of steady rain 
(solid line) reflects progressive washout, which was begun at 1 5, 30, 
and 50 km downwind. 
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Figure 26.- Dispersive decays of volumetric-average HCA concentration calculated 
as functions of downwind distance for Titan III SRM exhaust clouds for seven 
standard meteorological conditions at Cape Canaveral, Florida. 



Figure 27.- Dispersive decays of volumetric-average concentration calculated 

as functions of downwind drift time for Titan III SRM exhaust clouds for 
seven standard meteorological conditions at Cape Canaveral/ Florida. 
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Figure 28.- Calculated decays (solid line) of in-cloud concentration (verti- 
cal average, cloud center line) are shown for seven standard meteorologies at 
Cape Canaveral, Florida. Empirical fits of experimental measurements (dashed 
line) of peak (=1 sec) in-cloud HCA are shown for eight Titan III launches 
at Cape Canaveral. Note that various calculated and experimental HC?' decays 
are not canparable on a one-to-one basis, since respective input and actual 
launch meteorologies differed (from ref. 43) . 
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Figure 29.- Results of check on MDM-5(II) calculations of HCA(g) washout from 
SRM exhaust clouds. Check examines equivalency (circles versus straight 
line) of calculated HCit{g) retention in FFW Titan III exhaust cloud as 
determined by (a) ratio of calculated vertical HCA(g) column densities 
for respective rain and no-rain cases and (b) present simplified washout 
model. 


87 



Fractional retention of HCl(g) in SRM cloud, p(HC£)/p (HC£) 



Distance from cloud stabilization, km 

Figure 30.- Results of check on MDM-5(II) calculations of HCJl(g) washout from 
SRM exhaust clouds. Check examines equivalency (circles versus solid 
line) of calculated HCfi/(g) retention in SFW Titan III exhaust cloud, using 
same procedure as in figure 29. 
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